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NOMENCLATURE 

ac=  ^acxjacy"acz^  Acceleration  of  boat  C.G. 

aiJ  C0®fflc^ent  of  the  elastic  deflection  at 

point  1  due  to  unit  force  at  point  j. 

Co  ^an  Cbord  measured  normal  to  hinge  line  of  that  portion 

of  control  surface  and  the  area  that  lies  aft  of  hinge  line. 

60  'jlS^a^e  from  mass  ceoter  of  control  surface  to  hinge 

line  (see  table  of  subscripts  below). 

def orma tion  dlSplaCement  of  ^  dement  due  to  elastic 

h  H?r°nenh°f  perturbatlon  acceleration  from  equilibrium 
flight  position  of  an  element  of  boat  volume  dr. 

S  Acceleration  of  gravity. 

Ep  *  deformation.  °E  dlSplacement  of  dement  due  to  elastic 

Et  x-component  of  perturbation  acceleration  from  equilibrium 
flight  position  of  an  element  of  boat  volume  dt 

h  Depth  of  hydrofoil  -  positive  down. 

h  ~  hy,hz ^  Angular  momentum  vector  of  boat. 

V  (hxRdhyR'hzR)  Vector  angular  momentum  of  rotor  relative  to 

boat  axes. 

Orthogonal  unit  vectors  in  direction  of  x,y,z  axes  re¬ 
spectively. 

Reference  length,  usually  the  distance  between  two  convenient 
points  on  the  forward  and  aft  hydrofoil  systems.  For  the 
"  .  ,  °f  u^lf°rmity  it  is  recommended  that  the  longitudinal 
L  an?®  bet”een  the  forward-most  points  of  the  forward  and 

ofroo°t  cSr  ^  US6d;  diStanCg  b6tWeen  leadi^  ^ge 


HYDRONAUTI CS ,  Incorporated 


i 
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Distance  from  boat  C.  G.  to  control  stock 
(see  table  of  subscripts  below). 

Mass  of  boat. 


(see  Figure  2), 


m  Diagonal  mass  distribution  matrix. 

mc  Total  mass  of  control  surface  (see  table  of  subscripts  below). 
n  =  (Wn2)  Unlt  vector  normal  outward  to  surface  element. 


p  Pressure, 

p  Perturbation  in  P. 

p  Time  rate  of  change  of  perturbation  of  P. 

q  Perturbation  in  Q„ 

4  Time  rate  of  change  of  q. 

r  -  (x,y,z)  Position  vector  of  a  point  on  the  boat  relative  to 
its  C.G. 

r  Perturbation  in  R. 

r  Time  rate  of  change  of  r. 

s  Variable  of  Laplace  transform  (see  Equation  [58]). 
t  Time. 


u  Perturbation  in  U. 

u  Time  rate  of  change  of  u. 

v  Perturbation  in  V. 

v  Time  rate  of  change  of  v. 

w  Perturbation  in  W. 


w 


x.y,z 


Time  rate  of  change  of  w. 

The  boat  stability  axes.  A  right-hand  orthogonal  system 
of  moving  axes,  fixed  in  the  boat.  The  z-axis  is  directed 
toward  the  bottom  of  the  boat,  the  xz  plane  lies  in  the 
vertical  plane  of  symmetry  of  the  boat,  the  origin  o  is 
located  at  the  center  of  mass  of  the  boat  and  the  x-axis 
is  forward  and  coincides  with  the  direction  of  motion  when 
the  boat  is  in  equilibrium  flight. 
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WZa 


x  ,y  ,z 
o  o  o 


A=  A. 


B=  B.  . 
_  ij 


C=  C, 


A  system  of  axes  similar  to  the  x,y,  z  axes  except  that 
the  axes  lie  along  the  principal  axes  of  inertia  through 
the  mass  center  of  the  system. 

A  right-hand  orthogonal  system  of  axes  fixed  relative 

to  the  equilibrium  water  surface.  The  x  y  plane  is 

o  o 

fixed  parallel  to  the  equilibrium  plane  of  the  free  water 

surface,  ohe  z  direction  is  vertically  down  and  the  x 
u  o 

axis  lies  in  the  general  direction  of  the  initial  motion 
of  the  boat. 

Reference  area,  usually  projected  area  of  foil  system 
on  x-y  plane  of  body  axes  at  design  condition. 

Area  of  that  portion  of  control  surface  and  tab  area 
that  lies  aft  of  hinge  line. 

Hydrodynamic  added  mass  square  matrix  arising  from  ac¬ 
celeration  motion  in  the  elastic  degrees  of  freedom. 

Hydrodynamic  force  square  matrix  arising  from  velocity 
motion  in  the  elastic  degrees  of  freedom. 

Hydrodynamic  force  square  matrix  arising  from  torsion 
in  the  elastic  degrees  of  freedom. 


D=lt 


Differential  time  operator. 


Differential  time  operator. 


Du  Su  *  Typical  partial  derivative  operator  with  respect  to  u 


Moment  applied  to  control  by  all  power  sources  (see 
table  of  subscripts  below). 

Generalized  force  driving  the  n-th  normal  mode  (see 
Equation  [66]). 
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F 


G  = 


a 


Zw 


H 


H 


cu 


I 

I 


n 


x 


xa 


y 


za 


zx 


J 


K 


^yJ  Fz  )  Externally  applied  force  vector. 

(K,M,N)  Externally  applied  moment  vector. 

Typical  hydrodynamic  transfer  function  (see  Equation  [58]). 

Hydrodynamic  moment  about  control  stock  (see  table  of  sub¬ 
scripts  below). 

Typical  rate  of  change  of  control  hinge  moment  with  u<  the 

dH, 

partial  derivative  -  (see  table  of  subscripts  below);  (see 
note  below ) . 

Unit  diagonal  matrix. 

Effective  momt.it  of  inertia  of  the  control  system  relative 
to  control  stock  (see  table  of  subscripts  below). 

Generalized  inertia  in  n-th  mode  (see  Equation  [67]). 

Moment  of  inertia  of  boat  about  its  stability  x-axis. 

Moment  of  inertia  of  boat  about  its  principal  x  -axis 

a 

Moment  01  inertia  of  boat  about  its  stability  y-axis. 

Moment  of  inertia  of  boat  about  its  stability  z-axis. 

Moment  of  inertia  of  boat  about  its  principal  z  -axis 

a 

Product  of  inertia  of  boat  with  respect  to  its  zx-axes. 

Moment  applied  to  control  surface  by  manual  power  (see  table 
of  subscripts  below);  (see  note  below). 

Hydrodynamic  moment  component  about  the  stability  x-axis 
(positive  from  y  to  z)  -  rolling  moment. 
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1 


K<e 

2 

M 


M 


(e) 


N 


N 


(e) 


1 


P 

P 

ay 

P 

c 


P 


ex 


P 


rx 


P 

rz 


The  contribution  of  the  elastic  degrees  of  freedom  in  the 
-Lateral  modes  to  the  total  roll  moment. 


The  contribution  of  the  elastic  degrees  of  freedom  in  the 
lateral  antisymmetric  modes  to  the  total  roll  moment. 

Hydrodynamic  moment  component  about  the  stability  y-axis 
(positive  from  z  to  x)  -  pitching  moment. 


The  contribution  of  the  elastic  degrees  of  freedom  to  the 
total  M-moment. 

Hydrodynamic  moment  component  about  the  stability  z-axis 
(positive  from  x  to  y )  -  yawing  moment , 

The  contribution  of  the  elastic  degrees  of  freedom  in  the 
lateral  modes  to  the  total  yaw  moment. 

Component  of  o>  about  x-axis  (see  note  below). 

Product  of  inertia  term  for  ailerons  (see  Equation  [30]). 

Moment  applied  to  control  surface  by  power  sources  other 
than  manual  (see  table  of  subscripts  below). 

Product  of  inertia  of  elevator,  defined  in  Equation  [24]. 
Moment  of  inertia  term  for  rudder  (see  Equation  [27], 
Product  oi  inertia  term  for  rudder  (see  Equation  [27]. 


Q  Component  of  cx>  about -y-axis  (see  note  below). 

R  Component  of  co  about  z-axis  (see  note  below). 

T  Kinetic  energy. 

U  Elastic  strain  energy. 

U  Component  of  Uq  along  x-axis. 

— ► 

Uo=  (U,V,W)  Velocity  vector  of  boat  center  of  gravity. 
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V 

W 

w 

X 

X 


H 


X 


'Hi 


r(e  ) 


H 


JHi 
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Volume 

Time  rate  of  change  of  V 
Weight  of  boat 

Component  of  along  z-axis. 

s:oep:8itive  airection  °f  the  «**>^ 

Hydrodynamic  component  of  force  per  unit  volume  in  x-direction. 

x-dScUon?  C0mp0nent  °f  foroe  on  the  1-th  element  in  the 
Rate  or  change  of  X  hydrodynamic  force  with  control  deflec- 

^  ■■  ax 


z. 

1 


tion  6j  the  partial  derivative^ — - 

56  ‘ 

yH2xi0f-aSterfa°refo;ce!he  dlreCti0"  °P  utility 

st,p"t2™«Sn"otseei:sv?sn: of  freedom  in  the 

srx?r:^tSr^cette  positive  toipp  °f  ^mty 

Hydrodynamic  component  of  force  per  unit  volume  in  z-direction. 
SSrSloi?  C0mp°"ent  Df  rorCe  °"  the  1-th  element  in  the 
Net  Z-force  at  element  i. 


Net  Z-force  at  element  j . 


’R 

(e) 


The  contribution  of  the  rigid  body  motions  to  the  Z-force 

tntei:07t-1DUtl0n  °f  the  elastlc  degrees  of  freedom  to  the 
total  Z-iorce. 
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8 

c 


e 


e 

n 


(t 


) 


T) 

R 

0 


$ 

h 


p 

p 

p 


B 

c 


T 


I 

CD 

0) 

n 


4> 


Control  surface  rotation  angle  about  hinge  line  (see 

Figure  2  and  note  below):  (see  table  of  subscripts  be¬ 
low). 

Angle  between  the  principal  x^axis  and  the  stability 
x-axis  positive  when  the  x^axis  is  pointed  above  the 
x-axis  forward  of  the  C.G. 

Generalized  coordinate  giving  displacement  of  n-th 
normal  mode  (function  of  time  t). 

Perturbation  in  rudder  angle  (see  note  below). 

Time  rate  of  change  of 

Perturbation  in  elevator  angle  (see  note  below). 

Time  rate  of  change  of  t). 

Pitch  angle  (see  Section  1-3  and  Figure  l). 

Sweep  angle  of  control  surface  hinge  line. 


Perturbation  in  aileron  angle  (see  note  below). 

Time  rate  of  change  of 

Mass  density  of  water,  slugs  per  cu.  ft. 

Average  mass  density  of  an  element  of  boat  volume  dr. 


Normal  distance  to  mass  element  from  hinge 
trol  surface  (see  Figure  2),  (see  table  of 
below ) . 


line  of  con- 
subscripts 


Volume  of  boat;  dr  volume  element. 

Roll  angle  (see  Section  I~3  and  Figure  l). 

Angle  of  heading  (see  Section  I- 3  and  Figure  l). 

(Pj  QjR)  =  Angular  velocity  vector  of  boat  (see  note  below). 
Natural  frequency  of  n-th  mode  (Rad. /sec.). 


Perturbation  in  © 
Perturbation  in 
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MATRIX  NOTATION 


1' 

z  ' 

1 

(  } 
[  1 


Unit  row  matrix  which  operates  so  as  to  sum  all  the  terms 
of  the  column  matrix  on  its  right. 

Row  matrix  of  z. 


Column  matrix. 
Square  matrix. 


STATIC  DERIVATIVES 


X  ,X  ,X  ,X, 
u  v  w  la 


Y  ,Y  ,  Y  ,Y 
u  v  w  h 

2  ,Z  ,Z  .  Z, 
u.3  v3  vi3  h 

K  ,K  ,K  ,  K 
u"  V  w  h 

M  ,  M  ,M  j  M, 
u  v3  w  h 

u3  v3  w  h 


Rate  of  change  of  X  hydrodynamic  force  with  u,v,w 
or  h  at  steady  state  equilibrium  condition;  the 

partial  derivative  ,  ~  ,  M  or 


)v_ ^  Sw 


Sh 


ROTARY  DERIVATIVES 


W  WX* 


WVVY« 

wwv 

VVW^ 

WWh 

VWVN» 


Rate  of  change  of  X  hydrodynamic  force  with  p,q, 
r,D  or  <t>  at  steady  state  equilibrium  condition;  ’ 

the  partial  derivative;  -1^  ^  or  ^ 

'  bq3  dr3  be  5<t>  • 
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ADDED  MASS  DERIVATIVES 


x.. 

x.,x. 

u 

v  w 

Y.  , 

Y.,Y. 

u 

v  w 

z.. 

Z.,Z. 

u 

v  w 

x.. 

X.,X. 

P 

q  r 

Y., 

T.  ,  Y. 

P 

q  r 

z.. 

Z.,Z. 

P 

q  r 

Rate  of  change  of  X  hydrodynamic  force  with  u, 

ax  ax  ax 

v,  or  w  ,  ^  '  or  aw  ' 


ADDED  MOMENT 


K.  , 

K. 

,  K. 

u 

V 

w 

M., 

M.  . 

,  M . 

u 

V 

w 

N., 

N.  . 

>N. 

u 

V' 

w 

K., 

K. , 

,K. 

P 

q 

r 

M. , 

M. , 

M. 

P 

q 

r 

N.  , 

N., 

N. 

P 

q 

r 

OF  INERTIA  DERIVATIVES 

Rate  of  change  of  K  hydrodynamic  moment  with  u 

v  or  w  ;  the  partial  derivative  ,  —  or  ~ 

ou  dv  aw 


SUBSCRIPTS 


bar,  — 
i 

a 

e 

r 

n 


=  square  matrix 
=  initial  reference  value 
=  aileron 
=  elevator 
=  rudder 

=  the  n-th  of  a  series 
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SUPERSCRIPTS 


dot,  • 
prime,  1 
bar,  — 
arrow,  -*■ 


=  derivative  with  respect  to  time 

=  represents  non-dimensional  quantity  (see  Table  I) 
=  Laplace  transform;  also  mean  value 
=  vector 


Note  -  All  rotations  (control  surfaces,  angular  velocity  components, 
etc. )  and  moment  components  are  positive  according  to 
the  usual  right-handed  convention;  i.e.,  from  the  x- 
■axis  to  the  y-axis  to  the  z-axis  to  the  x-axis. 
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ABS TRACT 


The  hydrofoil  boat  in  a  seaway  is  a  very  compli¬ 
cated  system.  In  its  simplest  form  it  may  be  treated 
as  a  rigid  body.  However,  this  representation  is  not 
always  correct.  In  some  cases  it  may  be  necessary  to 
take  into  account  other  factors  such  as  the  gyroscopic 
effects  of  propellers,  the  dynamics  of  the  control  systems 
and  the  elastic  behavior  of  the  boat.  In  this  report 
the  general  equations  of  motion  of  the  hydrofoil  boat, 
treated  as  a  rigid  body,  are  first  derived  and  discussed. 
The  gyroscopic  effects  of  the  propellers,  the  dynamics 
of  the  control  systems  and  the  effects  of  the  elasticity 
of  the  boat  are  then  considered  and  related  to  the  equa¬ 
tions  of  motion  of  the  boat. 
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EQUATIONS  OF  MOTION  OF  A  HYDROFOIL  BOAT 
INTRODUCTION 

In  order  to  study  the  effect  of  various  parameters  on  the 
stability,  control  and  motion  in  a  seaway  of  hydrofoil  craft,  it 
is  necessary  to  solve  a  set  of  mathematical  equations  which  re¬ 
late  the  physical  properties  of  the  hydrofoil  boat  with  the  hy¬ 
drodynamic  forces  resulting  from  its  own  motion  and  that  of  the 
seaway.  The  fundamental  equations  used  in  such  analyses  appear 
in  various  sources  and  with  varying  degrees  of  generality  often 
without  derivation.  It  is  the  purpose  of  this  report  to  pre¬ 
sent,  in  one  place,  a  collection  of  the  basic  equations  of  motion 
of  the  hydrofoil  boat,  together  with  brief  derivations;  but  in 
sufficient  detail  to  afford  a  thorough  understanding  of  these 
equations.  No  attempt  is  made,  however,  to  present  here  methods 
of  estimating  the  hydrodynamic  forces. 

In  the  first  part  of  this  report  the  general  equations  of 
motion  of  the  boat,  treated  as  a  rigid  body  with  spinning  rotors, 
are  derived  in  terms  of  a  set  of  axes  fixed  in  the  boat.  These 
are  six  non-linear  scalar  equations  representing  the  components 
of  dynamic  equilibrium  in  each  of  the  six  degrees  of  freedom. 

Six  additional  equations  relating  the  motions  of  the  boat  in  body 
axes  to  the  orientation  and  motion  of  the  boat  with  respect  to 
fixed  axes  are  then  derived.  The  above  mentioned  set  of  twelve 
non-linear  equations  characterize  the  rigid  boat  with  spinning 
rotors  and  are  sufficient  to  determine  its  response  to  an  arbitrary 
set  of  time  dependent  forces  and  moments.  In  addition  the  equa¬ 
tions  characterizing  the  dynamic  response  of  each  of  the  major 
control  systems  are  derived. 
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The  second  part  of  this  report  discusses  the  applications 
of  the  foregoing  set  of  equations  to  the  determination  of  the 
stability,  control  and  response  to  a  seaway  of  a  hydrofoil  boat. 

Large  cross  coupling  terms,  in  these  equations,  between  the 
various  degrees  of  freedom  suggest  the  need  for  solving  all  the 
equations  of  motion  simultaneously  on  high  speed  digital  computers 
when  large  motions  such  as  those  obtained  in  a  heavy  or  sometimes 
even  mild  following  seas  are  of  interest.  Where  only  longitudinal 
motions  are  of  interest  considerable  simplification  may  be  achieved. 
Further  for  the  class  of  problems  in  which  the  boat  perturbations 
from  equilibrium  are  small  the  equations  of  motion  may  be  linear¬ 
ized  and  thus  lead  to  even  greater  computational  simplification 
with  acceptable  accuracy.  Such  linearization  also  makes  the  study 
of  automatic  controls  and  hydroelastic  effects  more  tractable. 

The  third  part  of  this  report  presents  the  linearized 
equations  of  motion  of  the  rigid  hydrofoil  boat  and  their  deriva¬ 
tion.  These  equations  are  also  derived  in  non-dimensional  form 
and  in  terms  of  hydrodynamic  transfer  functions.  The  latter 
set  of  equations  are  of  particular  importance  when  studies  of 
automatic  controls  and  effects  of  unsteadiness  on  boat  stability 
and  motions  are  of  interest. 

The  fourth  part  of  this  report  discusses  some  of  the  methods 
which  may  be  used  for  determining  the  effect  of  the  elastic  de¬ 
formation  of  the  boat  on  its  stability  and  control.  The  best 
known  of  these  are  the  method  of  normal  modes,  the  collocation 
approach  and  the  method  of  assumed  modes.  The  linearized  equations 
of  motion  of  the  boat  including  its  motion  in  the  elastic  degrees 
of  freedom  are  derived  by  the  first  two  methods.  A  set  of  these  three 
simultaneous  linear  equations  are  obtained  corresponding  to  the 
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original  three  rigid  body  degrees  of  freedom  in  the  longitudinal 
plane  (or  lateral  plane)  plus  n  degrees  of  freedom  in  the  elastic 
modes.  For  the  case  where  the  natural  frequencies  of  vibration 
of  the  structure  are  much  higher  than  the  frequencies  of  the 
rigid-body  motions  the  hydrodynamic  loading  may  be  considered 
slow  enough  so  that  the  structure  is  always  in  static  equilibrium. 

This  simplification  leads  to  the  quasisteady  equations  of  motion. 

The  steady  flight  equilibrium  equations  in  which  the  hydroelastic 
effects- are  taken  into  account  are  then  given  as  a  special  case 
of  the  former  set  of  equations. 

Much  of  the  material  has  been  adapted  from  standard  references. 
In  particular.  Reference  14  has  served  as  a  valuable  guide. 

_ DERIVATION  OF  EQUATION'S  OF  MOTION  OF  RIGID  BOAT 

1.  Rigid  Body  Equations 


In  this  Section  the  hydrofoil  boat  is  considered  as 
a  rigid  body  and  a  derivation  of  its  equations  of  motion  from 
first  principles  is  presented. 

From  Newton's  second  law,  the  force acting  through  the 
center  of  gravity  c  of  the  boat,  in  vector  notation  is 


dU 

"d  o 

F  =  m  — — 
dt 

where  m  is  the  mass  of  the  boat  and  its  velocity.  The  moment 
of  the  boat  1?  is  equal  to  the  time  rate  of  change  of  its  angular 
momentum  and  is  given  by 

7t  _  dti 
*  dt 

where  the  angular  momentum  "h  is,  in  terms  of  a  particle  of  mass 
dm  located  a  vector  distance  r  from  the  center  of  gravity, 

h  =  /¥  [  r  x  (  +  ax  r  )  ]  dm 


[1] 


[  2  J 
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where  the  integration  is  taken  over  the  whole  boat.  Since  "tJ 
is  the  same  for  all  particles  and  /  "r  d  m  =  0  the  above  equa¬ 
tion  becomes,  after  expanding  the  vector  triple  product, 

h  =  /  [  cjo  (r.~r)  -  ~r  (a>*r)  ]  dm 
Let  r  and  be  expressed  by 


r  =  xi  +  yj  +  zk 

oj  =  pT+q“j  +  R1? 


where  x,  y  and  z  are  the  scalar  components  of  ~r ,  P,  Q  and  R  are 
the  scalar  components  of  and  T,  1  and  Tc  are  unit  vectors  in 
the  direction  of  x,  y  and  z.  Substituting  these  expressions 
in  Equation  3  gives  for  the  scalar  components  of  7i,  in  terms 
of  the  moments  and  products  of  inertia  A,  B,  .  .  .  F, 


where 


h. 

x 

h 


y 

h 


z 


IP-I  Q  -  I  R 
x  xy  ^  zx 


-I  P  +  I  Q  -  I 

xy  y  yz 


I  P  -  I  Q  +  I 

zx  yz  z 


R 

R 


I 

X 

=  I 

(y2 

+ 

z2 

)  dm. 

I 

yz 

=  / 

yz 

dm. 

I 

y 

=  I 

(z2 

+ 

X2 

)  dm. 

I 

zx 

=  / 

zx 

dm. 

i 

Z 

=  / 

(x2 

+ 

y2 

)  dm. 

and  I 

xy 

=  / 

xy 

dm 

If  x,  y  and  z  are  taken  as  the  coordinates  of  dm  in  a 
non-rotating  frame  of  reference,  with  origin  at  the  center  of 
gravity  of  the  boat,  C,  then  it  is  clear  that  in  general  the 
moments  and  products  of  inertia  as  well  as  the  angular  velocity 
components  P,  Q  and  R  will  vary  with  time  as  the  boat  rotates. 
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This  is  an  unnecessary  complication  which  can  be  avoided  if  the 
coordinate  system  0  x  y  z  is  fixed  in  the  boat  and  allowed  to 
rotate  with  it.  Though  this  introduces  additional  terms  the 
resulting  equations  are  much  simpler  since  now  the  inertia 
terms  will  remain  constant.  Thus  Equation  [l]  becomes  in  terms 
of  the  velocity  components  of  Uq  (U,V,W)  and  angular  velocity  To 


F  =  m 


dt 


(u 


1  + 


V  j  +tf  k) 


=  m 


=  m 


(u  T  +  v"j  +  w  Tc)+  (u 


d  t 


L  +  W 


d  k, 
d  t ; 


6  U 

o  — ►  — ► 

-T—  +(D  X  U 
ot  o 


[5] 


since  dt-(DXi  ,  |-f  =  7  x  1  and  =  To  x  7  .  The  operator 

g 

Bt  has  the  definition  implied  by  Equation  [5].  Similarly,  Equa 
tion  [2]  becomes 


d  =  —  +  a)  X  h  . 

Combining  Equations  [5]  and  [6]  with  Equation  [3]  gives  for 
the  scalar  components  of  7  in  the  7,  7  and  7  directions 

Fx  =  m  [U  +  QW  -  R  V] 

Fy  =  m  [V+RU  -  PW] 

=  m  [W  +  P  V  -  Q  U] 


[7] 
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and  for  the  scalar  components  of  "5  In  the  ~i,  and  "Ic  directions 


K  =  ft  +  Q  h  -Rh 
x  ^  z  y 


M  =  h  +  R  h  -  P  h 

y  x  z 

N  =  h  +  Ph  -  Q  h  • 

z  y  x 


Equations  [7]  and  [8]  are  the  Euler  equations  of  motion  of  the 
hydrofoil  boat. 

2.  Effect  of  Spinning  Rotors 


In  the  derivation  of  the  Euler  equations  the  boat 
was  assumed  to  be  a  rigid  body.  When  there  are  spinning  rotors 
attached  to  the  boat  such  as  propellers  and  engine  rotors, 
their  effect  on  the  boat  motion  may  sometimes  be  significant 
enough  to  require  their  inclusion  in  the  equations  of  motion. 

In  this  case  it  becomes  necessary  to  modify  Equation  [4]  to  in¬ 
clude  the  effect  of  the  angular  momentum  of  the  rotors.  It 
can  be  shown  that  the  total  angular  momentum  of  the  boat  Is 
simply  the  vector  sum  of  that  given  for  the  boat  treated  as 
a  rigid  body  and  the  angular  momentum  of  each  rotor  relative 
to  the  body  axis.  The  latter  is  determined  from  Equation  [4] 
by  interpreting  the  moments  and  products  of  inertia  therein  as 
those  of  the  rotor  with  respect  to  axes  parallel  to  0  x  y  z  and 
origin  at  the  center  of  gravity  of  the  rotor.  Thus  if  the 
resultant  of  the  relative  angular  momenta  of  all  rotors  is  h 

R 


with  components  h  ,  h  , 

XK  yK 

(Equation  [4])  becomes 


h  then  the  total  angular  momentum 
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h  = 
x 


I  P 

X 


I  Q-I  R  +  h  _ 
xy  ^  zx  xR 


h  = 

y 


I  P  +  I 
xy  y 


Q-I  R  +  h  _ 
yz  yR 


h  =  -  I  P-I  Q+I  R  +  h 

z  zx  yz  *  z  zR 

and  the  moment  equations  (Equation  [8])become 


K  =  h 

+ 

Q  h 

-  R  h 

+ 

Q  h  - 

P  h 

X 

z 

y 

zR 

yR- 

M  =  h 

+ 

R  h 

-  P  h 

+ 

R  h  - 

P  h 

y 

X 

z 

xR 

zR 

N  =  h 

+ 

P  h 

-  Q  h 

+ 

P  h  - 

Q  h 

z 

y 

X 

yR 

*  xR 

[91 


[10] 


where  the  additional  terms  are  known  as  gyroscopic  couples. 
For  example  suppose  the  axis  of  a  rotor,  having  a  moment  of 
inertia  I  and  angular  velocity  fi,  is  in  the  direction  of  body 
axis  CK,  then  hR  =  i  I  n  and  the  components  of  the  gyroscopic 

couple  are  0,  I  R  and  -  I  fi  Q. 


3.  Motion  of  Boat  Relative  to  Fixed  Coordinate 

Solution  of  Equations  [7]  and  [8]  gives  the  linear 
velocity  components  U,  V,  W  and  angular  velocity  components  P, 

Q,  R  relative  to  the  0  x  y  z  axes  fixed  in  the  boat.  To  ob¬ 
tain  the  motion  of  the  boat  center  of  gravity  it  is  necessary 
to  express  the  linear  velocities  relative  to  a  fixed  coordinate 
system. 

We  therefore  define  a  fixed  orthogonal  coordinate  system 
°  Xo  yo  Zo  in  which  the  x0  Plane  is  fixed  parallel  to  the 

equilibrium  plane  of  the  free  water  surface  and  the  z  direction 

o 
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is  positive  downward.  The  orientation  of  the  boat  axes  (XjyjZ) 
relative  to  the  fixed  axes  (VVZ0)  is  shown  ln  Figure  1.  It 

is  assumed  that  at  first  the  two  reference  frames  are  parallel. 
Then  the  orientation  of  the  boat  is  determined  by  considering 
the  following  three  rotations,  in  the  order  indicated,  where  all 
rotations  are  in  the  positive  direction. 

i.  Rotate  the  XQ,yQ  axes  about  zq  through  the  angle  of 
yaw  ¥  to  the  position  (x^,  yj.  Then  the  direction  cosines  be¬ 
tween  (xQ,yo,zo)  and  (x^y^z  )  are  given  by 


X 

y 

z 

1 

1 

i 

x 

o 

cos  ¥ 

-  sin  ¥ 

0 

yo 

sin  ¥ 

cos  ¥ 

0 

z 

0 

0 

0 

1 

Thus  a  vector  x  ~i  +  y  ~j 

11  ii 

has  the  scalar  components 

x  =  x  cos  ¥  - 

o  1 

y  =  x  sin  ¥  + 
o  1 


+  z  k  in  the 

i  i 

y  sin  ¥ 

i 

y  cos  ¥ 

i 


x  , 

i 


y  . 

i 


z  system 

i 


in  the  (x0>yo>z0)  system.  Expressed  in  matrix  notation  these 
equations  may  be  written  as 
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cos  y  -  sin  I  o 

sin  ¥  cos  I  o 

0  0  1 

Thus  Equation  [lla]  represents  the  first  rotation. 

ii.  Rotate  the  (x^zj  axes  about  the  y  axis  through 

an  angle  of  pitch  0  to  (x^y^z^).  This  may  be  expressed  in 
matrix  notation  by 

cos  ©  0  sin  0  x 

2 

0  io  y  [nb] 

2 

-sin  0  0  cos  0  z 

2 

ill.  Rotate  the  (y  ,z  )  axes  about  the  x  axis  through 

2  2  2 

an  angle  of  roll  <J>  to  (x,y,  z)  the  actual  orientation  of  the  boat. 
This  may  be  expressed  by 

1  0  0 
0  cos  $  -sin  $ 

0  sin  $  cos  $ 

Combining  Equations  [lla],  [lib],  [lie]  gives  for  the  resultant 
of  all  three  rotations 


x 

y  [lie] 

z 


X 

0 

COS 

I 

-sin 

W 

0 

cos 

© 

0 

sin  © 

1 

0 

0 

X 

y0 

= 

sin 

I 

cos 

W 

0 

0 

1 

0 

0 

cos 

$ 

-sin$ 

y 

z 

0 

0 

0 

1 

-sin 

e 

0 

cos  0 

0 

sin 

cos$ 

z 
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By  carrying  out  the  operations  indicated  by  Equation  [12],  the 
direction  cosines  between  the  x,  y,  z  axes  and  the  x0*yo>z0  axis 
are  obtained.  These  are  tabulated  below. 


X 

y 

z  ■ 

X 

o 

cos  0  cos  Y 

-cos  $  sin  Y 
+  sin  9  sin  $  cos  Y 

sin  $  sin  Y 
- +  sin  0  cos  $  cos  Y 

yo 

cos  0  sin  Y 

cos  $  cos  Y 
+  sin  ©  sin  <E>  sin  Y 

-  sin  $  cos  Y 
+  sin  @  cos  $  sin  Y 

z 

o 

-  sin  0 

cos  0  sin  <£> 

cos  0  cos  $ 

Resolving  the  U,  V,  W  velocity  components  in  the  x0>yo>z0  direc¬ 
tions  gives  for  the  velocity  components  of  the  center  of  gravity 
of  the  boat  in  fixed  coordinates 


dx 

c 

dt 


dy^ 

dt 


dz 
_ c 

dt 


U  cos  0  cos  Y  +  V  (sin  0  sin  $  cos  Y  -  cos  $  sin  Y) 

+  W  (sin  sin  Y  +  sin  0  cos  $  cos  Y) 

U  cos  ©  sin  Y  +  V  (cos  $  cos  Y  +  sin  0  sin  $  sin  Y) 

+  W  (sin  0  cos  $  sin  Y  -  sin  0  cos  Y) 

~U  sin  0  +  V  cos  0  sin  +  W  cos  0  cos  $  [ 


To  obtain  x  y  z  from  integration  of  these  equations, 
o  o  o 

as  a  function  of  time,  in  the  most  general  case  would  clearly 
require  a  very  considerable  effort.  However,  for  most  problems 
Involving  the  motions  of  hydrofoil  boats  these  equations  may 
be  considerably  simplified  as  will  be  seen  later. 
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4.  Angular  Orientation  of  the  Hydrofoil  Boat 

The  angular  orientation  of  the  hydrofoil  boat  (e,0,Y) 
may  be  expressed  in  terms  of  the  angular  velocity  components 
(  P,  Q,  R )  by  expressing  the  angular  velocity  of  the  boat  in  terms 
of  the  angles  (0,O,T).  Thus  the  resultant  rotation  of  the  boat 
from  orientation  (@,0,?)  to  orientation  (#'+de,$+d<J>, Y+dY )  in 
time  dt,  using  the  same  order  and  axes  of  rotation  as  in  the 
previous  section:,  may  be  represented  by 

dft  =  (Y  +  dY)  Tc  +  (0+de)  7'  +  ($+do)  7 
°  i 

-  Y  7  -9  j  -  0  7 

°  x 

where  the  subscripts  indicate  the  axes  along  which  the  unit  vec¬ 
tors  7,  7  and  "k  are  taken  and  7  -«•  j  and  T'  —■  7  ad  dt  -*•  0 

i  i 


dfi  =  df  k  +  d©  j  +  dO  i 

O  i 

and 

7  =  §=y7o  +  ©7  +*7=p7+q7+r7  [  i4  ] 

Resolving  Y  k  and  ©  7  along  the  body  axes  x,  y,  z  by  the  use 

o  1 

of  Equations  [lib],  [lie]  and  [12]  gives 


P  =  &  -  Y  sin  0 

Q  =  6  cos  0  +  Y  cos  0  sin  $  [15] 

R  =  Y  cos  0  cos  0-0  sin  0 

Solving  Equation  15  for  0,  0  and  Y  gives 
0  =  Q  cos  0  -  R  sin  0 

=  P  +  Q  sin  0  tan  0  +  R  cos  0  tan  0 
Y  =  (Q  sin  ®  +  R  cos  0)  sec  0 


[16] 
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Integration  of  Equation  [ 1 6 ]  in  general  would  require  a 
complicated  numerical  procedure.  However,  for  many  problems  it 
will  be  shown  later  that  these  equations  may  be  considerably  sim¬ 
plified  without  serious  error. 

5-  Choice  of  Axes  Fixed  in  the  Boat 

The  equations  derived  in  the  preceding  sections  are 
valid  for  any  set  of  orthogonal  axes  fixed  in  the  boat  with  origin 
at.  the  center  of  gravity,  0.  However,  depending  on  the  problem  to 
be  solved  there  is  generally  one  orientation  that  is  to  be  pre¬ 
ferred  to  all  others.  Since  it  seems  fair  to  assume  that  all 
hydrofoil  boats  have  a  vertical  plane  of  symmetry,  it  is  generally 
desirable  to  make  this  plane  coincide  with  the  xz  plane  of  the  body 

axes.  This  makes  the  two  products  of  inertia  I  and  I  equal  to 

yz  xy 

zero,  thus  simplifying  Equation  [4],  This  leaves  the  directions 
of  the  x  and  z  axes  to  be  selected. 

Principal  Axes  -  If  the  x  and  z  axes  are  taken  to  co¬ 
incide  with  the  principal  axes  of  the  boat,  then  the  product  of  in¬ 
ertia  I  also  becomes  zero  further  simplifying  Equation  [4]. 

z_<  X 

These  axes  are  called  the  principal  axes  of  the  boat. 

Stability  Axes  -  If  the  x-axis  is  taken  in  the  direc¬ 
tion  of  horizontal  steady,  straight  flight  then  the  reference  ve¬ 
locity  components  V  and  W  are  zero,  thus  simplifying  the  equations 
of  motion.  The  z-axis  will  always  be  taken  as  positive  downward. 
These  axes  are  referred  to  as  stability  axes.  In  general,  these 
axes  do  not  coincide  with  the  principal  axes  and  will  obviously 
be  different  for  different  steady  flight  trim  conditions.  The  mo¬ 
ments  and  products  of  inertia  relative  to  these  axes  I  ,  I  and 

x  z 

I  will  therefore  depend  on  the  trim  conditions.  If  e  represents 
the  angle  between  the  principal  x  -axis  and  the  stability  x-axis 

3 
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and  the  principal  moments  of  inertia 
known  then 


I  and  I  for  the  boat  are 
xa  za 


I 

=  I  cos2 

e  + 

I  sin' 

X 

xa 

za 

I 

=  I  sin2 

€  + 

I  cos' 

z 

xa 

za 

I 

1 

H 

i— 1|(\ 

il 

I 

)  sin  2 

zx 

'  xa 

za 

[17] 


Body  Axes  -  When  the  x-axis  is  taken  parallel  to  some 
fixed  reference  line  in  the  boat,  such  as  the  line  from  which  the 
trim  angle  is  measured,  the  resulting  set  of  axes  are  called  body 
axes.  The  hydrodynamic  force  characteristics  of  the  boat  are 
usually  measured  with  respect  to  these  axes.  It  will  be  seen  later 
that  it  is  desirable  to  write  the  stability  equations  in  the 
stability  axes  system.  Therefore  it  will  in  general  be  necessary 
to  transform  all  the  force  and  moment  components  into  the  stability 
axes  system. 

6.  Equations  of  Motion  of  the  Control  Systems 

We  shall  assume  that  each  control  system  is  made  up  of 
a  control  surface  or  surfaces  (rudder,  elevator  or  ailerons)  as¬ 
sumed  to  be  rigid,  and  a  mechanical  linkage  of  rigid  elements  con¬ 
nected  to  a  power  source  (manual  and/or  hydraulic,  electrical, 
etc.).  Each  system  is  assumed  to  have  one  degree  of  freedom  for 
which  a  convenient  generalized  coordinate  is  the  control  surface 
angle.  The  motion  of  the  control  surface  will  be  determined  by 
the  moment  applied  by  the  power  source  E  ,  the  hydrodynamic  moment 

on  the  control  surface  H  ,  the  inertial  forces  of  the  system  and 

the  gravity  forces.  In  most  cases  the  gravity  forces  are  reasonably 
constant  and  not  important  so  that  they  will  not  be  included  in 
the  following  analysis. 
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The  forces  from  the  power  source  are  considered  to  be  made 

up  of  a  manual  control  force  J  and  a  power  control  force  P  which 

c 

result  in  an  applied  moment  on  the  control  surface  F  This  mo- 

c 

ment  is  related  to  Jc  and  Pc  by  the  linkage  system.  Thus  if  the 
displacements  due, to  Jc  and  Pc  are  6cJ  and  6cp  respectively  and 
this  results  in  a  rotation  6(6c)  of  the  control  surface  then  the 
work  done  by  Pc  and  Jq  neglecting  friction  are  given  as 


P  6(6)  =  P  6^  +  J  6 

c  c  c  Cp  c  Cj 


or 


d  „  d 

cP  +  J  cj 


F  =  P 

c  c  d6_  '  uc  d6 


tern  is 


c  ~’c 

The  inertial  force  of  a  mass  element  6  of  the  control  sys- 


:i8j 


6l"-5„?=-6m[5i|+?c+|^x7+2Sxl|  +  -x  [19] 


m 


where  the  operator  ^  is  defined  by  Equation  [5]  and  7  is  the  ac¬ 
celeration  of  6m  (see  Reference  l)  expressed  in  terms  of  r,  the  po¬ 
sition  vector  of  6^  in  the  boat  axes.  The  term 

—  6U  _ 

=  j—  +  at  x  U 
c  ot  o 

is  the  acceleration  of  the  center  of  gravity  of  the  boat.  Only 

that  component  of  6F'  that  acts  in  the  direction  in  which  the  mass 
element  is  free  to  move  relative  to  the  boat  axes  will  contribute 
to  the  motion  of  the  control  system.  The  work  done  by  6F 1  during 
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a  virtual  displacement  8  r  of  the  mass  element  is 

5P'  •  5  r  =  -6m6r  — —  +  [  "a^  +  X  r  +  a  x  (  uj  X?)  ]  •  6  ~r  [20] 

8t 

.  — ►  8  r  82  r 

since  8  r  is  colinear  with  — —  and  -  and  the  Coriolis  Force 

6t  6t2 

6m(2  cu  X  is  perpendicular  to  6  "r.  The  first  term  on  the  right 

of  Equation  [20]  is  the  inertial  work  on  8  due  to  its  acceleration 

m 

relative  to  body  axes  while  the  remaining  terms  are  due  to  the  ac¬ 
celeration  of  the  boat  itself.  Since  the  control  system  has  a 
single  degree  of  freedom,  the  magnitude  of  8  r  will  be  directly 
related  to  the  rotation  of  the  control  surface  8(5  )  and 

8F  •  (6  )  =  -  |  6  ~  + 

c  1  m  6(6o>  6t2 


— ►  o  cd  — ►  — ►  /—¥  — ► 

6  a  +  — —  X  r  +  a)  x  x  r 
m  c  8  t 


6(6C) 


where  8F  is  the  inertial  moment  at  the  control  surface  due  to  6  . 

m 

The  total  inertial  moment  is  obtained  by  summing  up  over  all  mass 
elements.  Thus  the  total  inertial  moment  at  the  control  surface 


where  I 


F  =  -I  - -  +  F. 

C  dt2  1 

dr  2 

-rrr — r  5  is  the  effective  moment  of  inertia  of 

d(6  )  m 


the  control  system  at  the  control  surface  and 
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\  e  d.  3?  r  5  CO  “►  — ►  /  —►  -4-  v  n 

A  5m  d[T)  aC  +  6~ t  Xr+o)X  (  a)  X  r  )  ] 

c 

X  5m  dfX)'  )  [ac  +  +  PR)Z  +  (^P-R)y  +  ( Q2+R2  )x]|i 

+  [a  +  (R+PQ )x+ (RQ-P )z- (P2+R2  )y ]  j 

y 

+  [ac  +.  (P+QR)y+(PR-Q)x-(p2+Q2)z]  "k  j  [23] 

Although  for  a  given  control  system  the  second  term  in  Equation 
[22]  may  be  computed.,  this  would  represent  a  considerable  effort 
except  for  very  simple  systems.  For  the  present  purpose,  we  will 
neglect  the  effect  of  the  control  linkages  in  the  computation  of 
this  term,  and  determine  the  effect  of  the  motion  of  the  boat  on 
the  control  surface  only.  The  control  surface  is  approximated  by 
a  thin  lamina. 

The  Elevator  System  -  Figure  2  represents  an  elevator 
free  to  rotate  about  the  swept  hinge  lines  shown.  It  is  assumed 
that  the  plane  of  the  elevator  is  parallel  to  the  xy  plane  of  the 
boat  axes  and  located  a  distance  &  aft  of  the  center  of  gravity 

and  a  distance  z  below.  A  small  mass  element  5  of  the  elevator 
^  m 

is  located  a  normal  distance  pg  from  the  hinge  line.  A  positive 
rotation  6g  of  the  elevator  is  considered  positive  for  the  trail- 

,  (J 

mg  edge  down.  Thus  d ,  >  =  pg  k  and  the  scalar  components  of  the 

'  °e  ' 

position  vector  of  Sm  is  given  by 

x  =  -  (ig  +  jy'j  sin  Ah  +  pg  cos  Ah),  y  =  y  and  z  =  zp  . 


F.  =  - 

1 
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Making  the  above  substitutions  In  Equation  [23]  gives  for  this 
term  alone 

Fi  =  ^acz  +  ^+!^)  y  ~  (pR-Q)  +  y'  sinA  +  p  cos  A,  ) 

J  B  n  e  h 

-  (P2  +  Q2)  zg  ]  pe  dm  . 

The  second  integral  is  zero  for  a  symmetrical  elevator.  Hence 

Fi  =  “  meee  [  acz“  (P2  +  Q2)ze  ]  "  Pex  (PR~^  [24] 

where  P0X  =  J (£^  +|y'|  sin  +  Pg  cos  A_h)  Pg  dm 

is  the  mass  of  both  elevators 
r  pedm 

eQ  )  Is  the  eccentricity  of  the  elevator  mass  center. 

J  e 

Equation  [24]  represents  the  inertial  force  on  the  elevator  due  to 
the  acceleration  of  the  boat.  This  effect  is  known  as  the  inertial 
coupling  of  the  elevator.  For  a  dynamically  balanced  elevator 
each  term  on  the  right  of  Equation  [24]  must  be  zero.  This  condi¬ 
tion  is  determined  by  the  following  consideration.  One  of  the 
required  conditions  is  that 

pex  =/^e  +iy'  slnA-h  +  Pe  cos  A'h  )  Pe 

=  ie  me  6e  +  Sln  Ah  [ pe  y'|  dm  +  cos  Ah  \  pe  dm  =  0 
With  the  other  required  condition  that  ee  =  0,  the  sweep  angle  of 
the  hinge  line  for  dynamic  balance  is  given  by 

Sln  Ah  j  pe  ly'|  dm  +  cos  Ah  !  P|  dm  “  0  [25] 

It  therefore  appears  that  by  "Statically  balancing  (e  =  0)  the 
elevator  about  the  hinge  line^  the  sweep  of  which  Is  given  by 
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Equation  [25],  it  is  possible  in  principle  to  dynamically  balance 
the  elevator.  In  practice,  however,  this  might  be  difficult  to 
achieve. 

Substituting  Equation  [24]  in  Equation  [22]  and  adding  the 
hinge  moment  applied  by  the  power  source  Fg  (see  Equation  [ 18 ] ) 

and  the  hydrodynamic  hinge  moment  H  ,  leads  to  the  equation  of 
motion  of  the  elevator  system. 
d26e 

Pe  “T  +  me  ee  [acz  '  <p2+«2)  ze]  +  Pex<PR-«)  =  Fe+He 

This  equation  does  not  include  the  effect  of  friction.  Where 
friction  is  a  significant  factor  the  effective  friction  moment 
at  the  hinge  must  be  estimated  or  measured  and  added  to  Equation 
[18]. 

The  Rudder  System  -  The  analysis  for  the  rudder  sys¬ 
tem  is  the  same  in  principle  as  for  the  elevator.  Figure  2  shows 
the  starboard  rudder  on  one  of  the  forward  struts  of  a  hydrofoil 
boat.  It  is  assumed  to  be  a  thin  lamina  parallel  to  the  xz  plane 
of  the  boat.  Its  hinge  line  is  taken  parallel  to  the  z-axis  lo¬ 
cated  a  distance  i  forward  and  y  to  starboard  of  the  boat  center 

of  gravity.  Therefore  we  have  for  an  element  of  mass  6  a  dis- 

m 

tance  aft  of  the  hinge  line 


X  =  £r  "  prJ  y  =  y 

Putting  these  conditions  into 
for  the  inertial  moment  on  the 


dr  ~t 

?  ’  d(iry  =  -  pr  J 

'  r 

Equa ti on  [23]  for  the  rudder  gives 
rudder  due  to  boat  acceleration 


Fj_  ~J  t  ac  y+  (R+PQ)  (^r -pr  h  (RQ-P)  Z  -  (P2+  R2 )  y^]  dm 
=  mr  er  [ac-y“  (?2+r2)  yrl  +  pr.x(R+pQ)  +  prz  (RQ-P) 


[26] 


[27] 
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where 

Prx=  1  ^ r~  Pr^  Pr  dm  is  a  moment  of  Inertia  term 

Prz=  f  PrZ  dm  ls  a  Product  of  Inertia  term. 

For  a  port  rudder  the  sign  of  y  is  changed.  For  a  rudder  whose 

hinge  line  is  a  distance  SL  aft  of  the  center  of  gravity  then 

?„  =  -  /  {l>  +  p  )  p  dm. 
rx  J  '  r  rr  r r 

The  equation  of  motion  of  the  single  rudder  system  is  an™ 
alagous  to  that  for  the  elevator  system 

d2  6r 

A  X“T  -  "V  e r  [acy-  (p2+R2)yr.]-PI,x(R+PQ)-PI,z(RQ-P)=Fr+Hr  [28] 


If  a  port  rudder  at  y  =  -  y^  is  rigidly  coupled  in  a  one  to  one 
ratio  to  the  starboard  rudder  then  the  term  containing  y  in 


Equation  [28]  cancels  out  for  the  coupled  system  and  the  equation 
of  motion  becomes 


d2  6 


r  p 

dtr 


2mreracy-  2Prx(R+PQ)-2Prz(R<l-P)  -  Fr  +  Hr 


+  Hr  [29] 

(starboard)  (port) 


where  1^  is  the  effective  moment  of  inertia  of  the  whole  rudder 
system  (including  both  rudders  in  the  latter  case). 


The  Aileron  System  -  The  aileron  system  differs  from 

the  coupled  rudder  system  in  that  one  control  surface  moves  up 

while  the  other  moves  down.  Although  these  surfaces  are  sometimes 

arranged  so  that  the  upgoing  aileron  moves  a  different  angle  than 

the  downgoing  one  it  will  be  assumed  here  that  they  move  through 

the  same  angle.  The  deflection  of  the  starboard  aileron  is  given 

Py  positive  downward  and  the  port  aileron  by  5;  positive  upward 

0. 
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The  ailerons  are  assumed  to  be  lamina  parallel  to  the  xy  plane 
with  hinge  lines  parallel  to  the  y-axis  which  are  located  a  dis¬ 
tance  i  forward  and  below  the  boat  center  of  gravity.  There¬ 
fore  we  have  for  an  element  of  mass  6^  a  distance  a  aft  of  the 
hinge  line 


x  =  Vv^’  Z=V  dfvr 


Putting  these  conditions  into 
face  and  adding  gives 


k, 


d  r 


a  (starboard  )>  d(6  ) 

a  ' 


Equation  [23]  for  each 


Pa  k(port) 
control  sur- 


~  2  f  (P  +  OR)  |y|  p  dm 

-L  1  a 

since  the  terms  in  a^,  x  and  z  cancel  out.  The  equation  of  mo 
tion  then  becomes 


I 

a 


d2  ( 6  ) 
a ' 

dt2 


+  2  P 


ay 


(P  +  QR)  = 


^a(SB)  ^a(Port) 


+  F 


[30] 


where  Ig  is  the  effective  moment  of  inertia  of  the  whole  .aileron 

i 

system  (including  both  surfaces) 

Fa  is  the  generalized  aileron  control  force  (see  Equation  [18]) 

Pay  =  1  |yl  Pa  dm  • 


7-  External  Forces 

The  forces  and  moments  on  the  left  hand  side  of  Equa¬ 
tions  [7]  and  [8]  are  the  external  actions  on  the  boat.  These  are 
in  general  of  three  kinds,  gravitational,  aerodynamic  and  hydro- 
dynamic.  The  gravitational  forces  are  determined  by  the  orienta¬ 
tion  of  the  boat.  Unless  the  boat  has  an  aerodynamic  propulsion 
system  the  aerodynamic  forces  will  usually  be  small  except  for  high 
speed  boats  in  which  case  wind  drag  is  generally  an  important  fac¬ 
tor.  The  hydrodynamic  forces  are  determined  by  the  elevation, 
orientation,  configuration  and  motion  of  the  boat,  z  ,  (e,  <t>,  l), 
(V  V  6a '  ‘  ^U'-  Vj  Wk'  (p>  Q-'  R)  and 
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the  characteristics  of  the  propulsive  system. 

Prom  Figure  3  the  components  of  weight  in  the  directions 
of  the  axes  are  readily  seen  to  be 
X.t  =  -  mg  sin  0 

fcj 

Y  =  mg  cos  0  sin  0 

o 

Z  =  mg  cos  0  cos  <t> 

O 

Thus  the  external  forces  in  Equation  [7l  may  be  written 
F  -  X  -  mg  sin  0 

A 

F  =  Y  +  mg  cos  ©  sin  $>  [31] 

F^  =  Z  +  mg  cos  0  cos  i> 

where  X,  Y  and  Z  are  aerodynamic  and  hydrodynamic  forces  (includ¬ 
ing  propulsive  forces).  Since  gravity  does  not  contribute  to 
the  moments  about  the  mass  center,  then  (L,  M,  N)  are  entirely 
aerodynamic  and  hydrodynamic. 

II.  APPLICATIONS  OF  EQUATIONS  OF  MOTION 
1.  Summary  of  Equations 

The  kinematic  and  dynamic  equations  of  motion  of  the 
hydrofoil  boat  and  control  systems,  derived  in  the  foregoing  sec¬ 
tions,  are  collected  below  for  ready  reference.  The  xz  plane  of 
the  coordinate  system  is  taken  in  the  plane  of  symmetry  of  the 
boat  so  that  Iyz  =  Ixy  =  0.  The  boat  is  assumed  to  be  a  rigid 

body  which  may  have  any  number  of  rigid  spinning  constant  speed 
rotors  attached  to  it. 

X  -  mg  sin  0  =  m  (U  +  QW  -  RV)  (a) 

Y  +  mg  cos  0  sin  T  =  m  (V  +  RU  -  PW )  (b)  [32] 

Z  +  mg  cos  0  cos  T  =  m  (W  +  PV  -  QU)  (c) 
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K  =  Ix  ^zx  -  W^zR  -  R1Yr  <a> 

m  =  Iy  4+RP(ix-iz)  +  izx(p2-R2)  +  RhxR-RhzR  (R)  1331 

N  =  Y  R-^x  +  Yx  ^R+PhyR  -  *XR  <°) 


P  =  $  -  ¥  sin.  6  (a  ) 
Q  =  0  cos  i>  +  ¥  cos  8  sin  $  (b) 
R  =  f  cos  0  cos  3>  -  0  sin  <5  (c) 
0  =  Q  cos  $  -  R  sin  cl  (d) 
<j>  =  P  +■  Q  sin  $  tan  0  4-  R  cos  0  tan  0  (e) 
¥  =  (Q  sin  $  4  R  cos  <J>)  sec  0  (f) 


[34] 
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Plications  of  Equations 


Equations  [34]  clearly  represent  only  three  independent 

equations .  Additional  equations  must  be  added  to  Equation  [36] 

when  more  than  one  control  system  of  a  given  type  is  used.  For 

example  if  the  boat  has  independent  fore  and  aft  elevators,  there 

would  be  a  separate  equation  for  each  elevator  system.  Thus,  in 

general,  for  n  independent  control  systems  Equations  [36]  consist 

of  n  independent  equations.  Since  the  aerodynamic  and  hydrodynamic 

forces  and  moments  (X,  Y,  Z),  (K,  M,  N)  and  H  depend  only  on  the 

n 

sea  state  and  the  elevation,  orientation,  configuration  and  motion 
variables  zq,  (6,  -5)  6^,  (U,  V,  W)  (P,  Q,  R)  they  do  not  repre¬ 
sent  independent  unknowns.  Also  Equation  [35a,  b]  may  be  dropped 
since  the  remaining  equations  do  not  depend  on  x  and  y  .  Thus 
Equations  [3?! -.36]  may  be  reduced  in  general  to  10  +  n  independent 
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equations  in  10  +  2n  unknowns.  The  unknowns  are  (U,  V,  W), 

(P,  Q,  R),  (8j  $,  W),  z  ,  6  and  F  where  5^  and  F  are  the 

°  cn  cn  cn  n 

angle  and  control  force  on  the  n-th  control  system.  It  is  there¬ 
fore  clear  that  n  of  these  variables  must  be  specified  before 
these  equations  can  be  solved  and  the'  motion  determined.  The 
variables  selected  and  the  manner  in  which  they  are  specified  are 
determined  by  the  problem  to  be  solved. 

Stability  Problems  -  Controls  Fixed  -  In  these  prob¬ 
lems  the  boat  is  initially  in  steady,  level,  flight  in  a  calm 
sea  and  it  is  desired  to  determine  the  motion  caused  by  a  dis¬ 
turbance  of  very  short  duration.  The  controls  are  assumed  locked 
in  the  position  for  steady,  level  flight.  Since  the  control 

angles  5  are  fixed  and  known  Equations  [36]  are  not  needed  to 
cn 

determine  the  boat  motion.  However,  after  the  boat  motion  is  de¬ 
termined  these  equations  may  be  used.  If  desired,  to  obtain  the 

control  forces  F  necessary  to  hold  the  control  surfaces  in  the 
cn 

fixed  condition.  Since  the  aerodynamic  and  hydrodynamic  forces 
do  not  depend  on  xq  and  y  ,  Equations  [35  a,b]  may  also  be  dropped 

The  result  is  a  set  of  ten  equations  in  ten  unknowns.  These  equa¬ 
tions  are  non-linear  and  difficult  to  solve.  For  many  practical 
problems  it  is  sufficient  to  linearize  the  equations  by  dealing 
with  small  perturbations  from  equilibrium  flight.  In  this  case 
the  resulting  equations  are  a  set  of  homogeneous  linear  differen¬ 
tial  equations  with  constant  coefficients  which  are  readily  solved 

Stability  Problems  -  Controls  Free  -  These  problems 

are  similar  to  those  discussed  above  except  that  in  this  case  the 

control  angles  are  variable  and  Equations  [36]  with  the  control 

forces,  F  ,  set  equal  to  zero,  must  be  added  to  the  above  system 
cn 
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of  equations.  If  the  equations  are  linearized  they  again  become 
homogeneous.,  and  are  treated  in  the  same  manner  as  in  the  fixed 
control  case. 

Stability  Problems  -  Automatic  Control  -  In  these 

problems  the  controls  are  neither  fixed  nor  free.  Even  though 

the  pilot  input  Pc  is  zero  the  generalized  control  forces  are 

n 

not.  They  are  in  general  determined  by  a  feedback  of  the  boat 
response  into  the  input  of  the  control  system.  All  10  +  n  equa¬ 
tions  are  in  general  involved,  and  the  dynamics  of  the  closed 
loop  must  be  introduced. 

Response  to  Controls  -  In  these  problems  it  is  desired 
to  determine  the  effectiveness  of  the  boat's  controls.  This  is 
normally  done  by  specifying  a  particular  variation  of  the  con¬ 
trol  surface  angles  6c  with  time,  e.g.  a  step-function  input  of 

n 

rudder  angle.  The  equations  of  motion  of  the  boat  then  become 
inhomogeneous  equations  for  (U,  V,  W),  (P,  Q,  R),  (© ,  <&,  y)  and 
zq.  The  control-systems  equations  are  not  needed  but  may  be  used 
if  desired  to  compute  the  variations  in  the  control  forces  F 

Ln 

required  to  produce  the  maneuver. 

Response  to  Seaway  -  The  response  of  a  boat  to  a  sea¬ 
way  is  treated  with  the  same  equations  as  are  used  for  stability. 
The  effect  of  the  seaway  introduces  new  forcing  functions  on  the 
lifting  surfaces  and  the  struts.  The  equations  of  motion  then 
become  inhomogeneous  and  the  mathematical  treatment  is  similar 
to  that  used 'for  the  response  to  controls.  For  seas  in  which 
the  motions  are  large  it  is  necessary  to  solve  the  equations 
without  the  simplification  of  linearization.  This  is  in  general 
a  formidable  task  requiring  a  high  speed  digital  computer  and 
has  heretofore  been  carried  out  in  a  few  cases  and  only  for  simple 
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foil  systems  performing  longitudinal  motions  in  regular  head  and 
following  seas.  It  has  been  found  that  when  the  motions  are  not 
very  large,  linearization  of  the  equations  of  motion  give  a 
reasonably  good  approximation  of  the  boat  motions.  It  has  been 
customary  to  use  the  linearized  equations  in  the  study  of  auto¬ 
matic  control  systems  in  relation  to  the  seaway. 

— — — r s e Isms  -  This  class  of  problems  arises  when 

specified  variables  are  selected  from  those  that  are  normally 

treated  as  dependent,  (U,  V,  W)  (P,  Q,  r),  (g,  O,  Y)  and  z  .  For 

o 

example  it  may  be  required  to  find  the  control  forces  required 
to  produce  a  certain  steady  level  turn.  In  this  case  one  might 
specify  the  roll  angle  3>,  the  turning  radius  R  (=  U  /Y),  the 

elevation  of  the  center  of  gravity  Zq(=  0)  and  n  -  3  control  sur¬ 
face  angles.  Since  the  motion  Is  assumed  steady  this  would  lead 
to  a  set  of  algebraic  equations.  If,  on  the  other  hand,  Y  and  i 
are  specified  as  functions  of  time  a  set  of  differential  equations 
are  obtained.  A  decided  advantage  of  this  method  is  its  ability, 
in  many  cases,  to  deal  with  non-linear  equations. 

III.  THE  LINEARIZED  EQUATIONS  OF  MOTION 

1.  Derivation  -  The  Dimensional  Equations 

The  solutions  to  the  equations  of  motion  when  the  mo¬ 
tion  of  the  boat  is  limited  to  small  infinitesmal  disturbances 
of  a  symmetrical  boat  from  steady  equilibrium  flight  in  a  calm 
sea  has  been  extensively  studied  by  many  investigators.  Although 
they  apply  strictly  to  infinitesmal  disturbances  only,  the  results 
obtained  have  been  found  to  apply  with  fair  accuracy  to  finite 
disturbances  which,  are  sufficiently  large  to  make  the  conclu¬ 
sions  drawn  of  great  practical  interest. 


HYDRONAUTICS,  Incorporated 


In  the  following  discussion  we  shall  carry  through,  the 
details  of  the  linearization  for  Equation  [32a]  and  [33a]  only. 
The  other  motions  equations  are  treated  in  an  analogous  manner. 
As  stated  earlier  the  aerodynamic  and  hydrodynamic  forces  and 
moments  depend  on  the  elevation,  orientation, configuration  and 
motion  variables.  Thus  in  Equation  [32a]  and  [33a] 

x  =  x  (U,V,W,U,V,W,P,Q,R,  P,Q,R,0,(D,T.6  ,  Be  .&  ,z  ) 

'  an  en'  rn  0 

K  =  K  (U,V,W,U,  V,W,P,Q,R,P,Q,R,0,O>,T  ,6.  .6  ,5  ,z  ) 

'  V  en  rn  0 

The  boat  is  assumed  to  be  In  some  equilibrium  reference  flight 

condition.  The  initial  reference  values  of  all  the  variables 

are  denoted  by  a  subscript  i,  and  the  small  perturbations  are 

indicated  as  follows:  changes  in  U  V  W  etc.  are  indicated  by 

small  letters,  i.e.  U=u  +  u,  P=p  +  p,e  =  0.  +  0,  z  =  z  +  h 

x  1  1  0  .1 

etc.  Changes  in  the  control -surface  angles  are  denoted  by 

(C,  r),  C  )  so  that 

8  —  61  +  P 

o 


,  6  ,  z  ) 

n 

[38] 

,5  ,  z  ) 

rn  0 

[39] 

8t  =  5  4-  Ti 

en  eni  n 


6=6  + 
rn  rni 

For  steady  equilibrium  the  reference  values  of  the  variables  are 
taken  as  constant.  Making  the  appropriate  substitutions  in  Equa¬ 
tions  [32a]  and  [33a]  yields 

x  -  mg  sin  (0^+0)  =  m  [u+(q_^+q)  (w^+w  )-(r  ,+r  )  (v^+v)  ] 

K‘  V'  V  +  <V<1>(ri+I’)(lz-Iy)-I2X(P1+p)(q1+q)+(q.+q)hzn 

-  (r,  +r  )h 

i  yn 


[4o] 


[42] 
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where  X  and  K  are  given  by 


+  U,  V.  +  V,  - 

1 

-  5 

r  . 

©  ni 

+  ^ 

z±  +  h) 

[43] 

+  U,  V.  +  V,  - 

1 

- ,  6 

rni 

+  cn. 

z±  +  h) 

[44] 

The  Taylor  expansion  of  Equations  [43]  and  [44]  Is 

(uD  +vD  H - +hD  ) 

X  =,  [1+(uDu+vDv+— +hDh)  +  - ^ - |y - —  +— ]  X(u.  ,  v±,  z±  ) 


(uD  +vD  + - +hD  ) 

K  =  [  1+(uDu+vDv+— -+hDh)  +  — — - - - 1- ]  K(u±J  v1,  — -j  z1  ) 


Since  we  are  dealing  with  small  perturbations  frdm  the  reference 
flight  condition,  the  second  and  higher  order  terns  may  be  omitted 
leaving  only  the  first  two  terms  of  each  of  the  above  series. 
Writing  the  above  operators  D  in  their  partial  derivative  notation 


X  (u. ,v.  -- 
i  i 


V  ,  dX 

■-  zi>  +  u  si 


+  v 


ax 

dv 


+ - +  h 


ax 

ah 


45] 


K  =  K  ( u . ,  v . 

l  i 


%  ,  ax 

Z  )  +  U  y— 
l  du 


+  V 


dK 

dv 


+ - +  h 


8K 

dh 


[46] 


where  the  subscript  i  denotes -that  the  partial  derivative  is  taken 
at  the  initial  equilibrium  condition.  The  use  of  only  first  order 
terms  implies  the  assumption  that  the  forces  and  moments  vary  lin¬ 
early  with  the  disturbance  variables  for  small  enough  disturbances. 
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Also  in  its  present  form  it  is  assumed  that  the  derivatives  de¬ 
pend  only  on  the  instantaneous  values  of  the  disturbance  velocities, 
control  angles,  etc.  and  their  derivatives.  Despite  the  success 
of  this  last  assumption  in  many  types  of  problems  it  is  in  general 
not  sound.  It  may  lead  to  serious  errors  in  cases  where  the 
hydrodynamic  forces  change  very  rapidly,  as  in  flight  in  rough 
seas  at  high  speeds  or  when  a  control  is  very  rapidly  displaced. 

An  alternative  set  of  small  disturbance  equations  will  be  given 
later  which  is  not  subject  to  these  limitations. 

On  substituting  Equations [ 45 ]  and  [46]  in  Equations  [4l] 
and  [42]  and  expanding  sin  (0.  +  e)  one  obtains 


v  ,  dx 

X.  +  U  'c — 
l  du 


+  v 


dx 

dv 


+  ---+  h 


dX 

dh 


-  mg  (sin  6.  +  Q  cos  0.  ) 
1  1  ' 


-  m  [u  +  (q^+q)  (w^+w)  -  (r_^+r)  (v.  +v)] 


[47] 


.7-  ,  d  K 
K .  +u 

i  du 


+v 


dK 


+ - hh 


dK 


Sh[-  IxP-Izxi'+(cli+q)(r-i+r)(lz-Iy) 


(Pj+p)  (q,+q)  +  (q,+q)  hzR-  (r,+r)  h 


ZX  '  1 


yR 


[48; 


Since  X. 

i 

hand  side 
cancelled 


for 


dx 

du 


i 


and  are  equal  to  the  unperturbed  terms  on  the  right 

of  Equations  [47]  and  [48]  respectively  they  may  be 
out.  Retaining  only  first  order  terms  and  writing  X 

u 

etc.  gives 


(X.-m)  u+X.v+X.w+X  u  +  (X  +mr. )v+( 
u  vwu  vi/v 


■\r 


A 

W 


+  (Xq-mw±  )q+(Xr+mv1  )r+(X0-mg  cos  0  1  Je+X^+X^+X^+X^fi 

+  X^+X  e+X  n  +X.C  +  X,  h  =  0 

C  t,  4  4  4  h 


[  49] 
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and 

V+Kvi+V+Kuu+Kvv+K»"+(K{,-Ix)p+V + 

+  [  Vri  ( Vb )  ■hzR+IzxPi  1  q+[  VWb  >+hyR ] V+’Vt+ 

i+K^q+K^  k+K^  |+K^rj+K^+K^h  =  0 


The  remaining  equations  in  the  set  of  Equations  [32]  -  [36]  may 
be  determined  in  an  analagous  manner.  These  equations  give  the 
perturbation  motions  from  any  equilibrium  flight  condition. 

For  the  very  important  case  of  straight,  level,  symmetrical 
flight  these  equations  may  be  considerably  simplified.  For  a 
truly  symmetrical  boat  it  is  clear  that  the  side  force  Y,  rolling 
moment  K,  yawing  moment  N,  and  aileron  and  rudder  hinge  moments 
and  will  be  zero.  Thus  the  derivatives  of  the  asymmetric 

forces  and  moments  with  respect  to  the  symmetric  variables  u,  w, 
cb  q,  0,  h  are  all  zero.  Furthermore  v.  =  p.  =  q.  =  r.  =  .  =0. 

t 

We  will  use  stability  axes  so  that  w^  =  0  and  u.  is  then  the  re¬ 
ference  boat  speed  U  .  Since  we  are  considering  level  flight 
e±  =  °-  In  addition  the  following  approximations  will  be  made. 


1.  We  may  neglect  the  derivatives  of  the  symmetric 
forces  and  moments  X,  Z,  M,  Hp  with  respect  to  the  asymmetric 

variables  v,  p,  r,  4>,  £. 

2.  We  may  neglect  all  the  acceleration  derivatives 


except  the  following  Z.,  Z^,  M.,  He.,  H^,  Y.,  Y.,  Y.,  Ky, 

h'  KrJ  NoJ  Kr’  Ha->  pr ■  ’  Hn.j  IC  •  In  practice  it  will  turn 

F  P  p  v  h 

out  that  most  of  the  above  derivatives  will  often  be  negligible. 


[50] 
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3.  We  may  neglect  H^,  Hg^  and  Hr^  and  the  angular 
momentum  of  the  rotors. 


When  the  above  simplifications  are  made  the  linearization  of  Equa¬ 
tions  [32]  -  [36]  for  calm  sea  operation  results  in  Equations 
[51]  and  [52]  (page  32).  The  effects  of  a  seaway  are  taken  into 
account  by  simply  adding  the  hydrodynamic  forces  and  moments  im¬ 
posed  by  the  sea  to  the  right  hand  side  of  these  equations.  In 
general  these  would  be  a  function  of  time  and  the  perturbation 
variables  of  the  boat.  However,  useful  results  have  generally 
been  obtained  by  neglecting  the  dependence  of  the  wave  forces  on 
the  boat  perturbations  as  long  as  the  boat  motions  are  not  too 
vio.ent.  In  this  case  the  wave  forces  and  moments  are  computed 
as  if  the  boat  were  moving  in  steady  level  flight  and  appropri¬ 
ately  added  to  the  right  hand  sides  of  Equations  [51]  and  [52]. 

Due  to  the  simplifying  assumptions  that  were  made  in  the 
previous  section.  Equations  [  51 ]  and  [52]  have  become  uncoupled. 
Equations  [51]  are  a  function  only  of  u,  w,  q,  0,  h,  x,  and  r), 
the  longitudinal  variables.  Solution  of  this  set  of  equations 
gives  the  perturbed  longitudinal  motions  of  the  boat.  Equations 
[52]  contain  only  the  variables  v,  p,  r,  <t>,  Y,  y,  i  and  So¬ 
lutions  of  these  equations  gives  the  perturbed  lateral  motions 
of  the  boat. 


In  the  computation  of  longitudinal  motions  it  is  often 


convenient  to  eliminate  w  from  Equations  [51a, b,c].  Substituting 
Equation;  [ 51 f ]  into  Equations  [51a,b,c]  gives 


(mD  Xu)u  ^XwI>fXh^h_^XqD4'XwUo4'X0"mg^  0=0 

-ZuU+[  (m-z.  )D*-ZD-Zh]h-[  Z.  D2+  (Zq+Z.  Uo  18  -Zq 

-MuU-(MliD2+MuIM-Mh)h+[  (ly-W4)D2-(Mq+MiiUo)D-(MwUo+M0)]8 


(a) 

3  =  0  (b) 


-(M  D+M  )n  =  0 

1  3 


(c)  [531 


longitudinal  equati ons 
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_ Dimens ionless  Equations 

Equations  [51a,b,c]  and  [52a,b.,c]  are  non-dimension- 
alized  by  dividing  the  force  equations  by  JrpAU^2  and  the  moment 

equations  by  ipAiu/.  The  quantity  A  is  a  reference  area  usually 

taken  as  the  projected  area  of  the  foil  system  on  the  xy  plane 
of  the  body  axes  and  t  is  a  reference  length  usually  the  distance 
between  two  convenient  points  on  the  forward  and  aft  hydrofoil 
systems.  For  the  sake  of  uniformity  it  is  suggested  that  the 
longitudinal  distance  between  forward-most  points  of  the  forward 
and  aft  foil  systems  be  used. 

Linear  velocities  are  non-dimensionalized  by  dividing  by 
Uq  and  angular  velocities  by  U/J .  The  mass  terms  are  divided 

by  *PAJ.  The  control  surface  equations  are  non-dimensionalized 

by  dividing  the  hinge  moment  H  and  control  force  F  by  4-pA  “  U  2 

°  c  “  C  C  O 

where  Aq  is  the  area  of  that  portion  of  the  control  surface  and 

tab  area,  that  lies  aft  of  the  hinge  line  and  C  is  the  mean  chord 

oi  tne  Same  portion  of  the  control  surface  and  tab  area.  The  non- 
dimensional  quantities  are  indicated  by  the  prime  symbol.  Thus, 
for  example 

X 

v  i  ,,  i  n  u 


K  ' 


P 


P' 


K 

-  ■  P  .sl 

!pi2AU  Uo 
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He  '  u' 
u 


He 
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ipFTl T 
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U 
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U  1  =  ~  =  1  ;  H  •  1  ?  '  - 
o  U  ’  r  £  ^  ~ 


i-g-  =  — gg. 

?pUo2  A 


=  C 


L. 


H 


lL 


ipc  A  U  i 
r  e  e  o 


D'  £ 


A  summary  is  given  in  Table  I. 


DIMENSIONAL 
QUANTITY  ’ 


X  Y  Z 
K  M  N 


TABLE  I. 

THE  NON-DIMENSIONAL  SYSTEM 
DIVISOR 


IpU  2A 
o 

IpU  2Ai 

O 


NON-DIMENSIONAL 
_ QUANTITY 


X'  Y>  Z' 
K'  M'  N' 


H  H  H 

a  e  r 

IpU  2A  6 
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H  ’  H  '  H  ' 
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F  F  F 
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IpU  2 A  c" 
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Equations  [51]  and  [52]  may  be  written  in  dimensionless 

form  simply  by  placing  primes  in  the  manner  above  indicated  on 

all  the  terms  and  making  the  substitutions  U  '  =  1  and  m'g!  =  CT  . 

o  ijo 

The  resulting  dimensionless  equations  are  sometimes  written  With¬ 
out  the  prime  symbols  for  the  sake  of  simplicity  in  notation. 

By  making  the  above  substitutions  in  Equations  [51],  [52],  [53] 
the  dimensionless  Equations  [54],  [55],  [56]  are  obtained  (see 
page  35). 


3 •  Hydrodynamic  Transfer  Functions 


In  the  previous  sections  the  hydrodynamic  derivatives 
were  written  in  quasistatic  form.  For  example  the  change  in  the 
dimensionless  Z' -force  due  to  a  perturbation  w'  was  written  as 


Z  '  =  (Z  '  +  Z.  '  D1  )  w'  [57] 

This  is  correct  as  long  as  w1  is  a  slowly  changing  function  of 
time.  However,  when  rapid  motions  are  involved  the  above  equation 
does  not  correctly  represent  the  force.  In  general  however  the 
relation  between  .Z  '  and  w'  may  be  assumed  to  be  representable 

by  a  linear  differential  or  integral  equation  of  more  complex 
form.  The  hydrodynamic  transfer  function  relating  w'  and  Z  '  is 


defined  as  the  ratio  of  the  Laplace  Transforms  of  Z  '  and  w 

1 


Gz»(s>  “ 


Z'(s) 

1 

w‘(s) 


:  58] 


-  00  -st1 

where  Z'(s)  =  /  e  Z'(t')  dt'  ;  etc.  (See  Reference  17) 

1  t'=0  1 

The  transfer  function  is  unique  if  w' (t)  =  Z  ' (t)  =  0  for  t'<  0. 

This  requires  that  we  deal  with  motions  for  which  the  initial  con¬ 


ditions  are  those  of  the  reference  flight  conditions.  Because  of 
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the  assumed  linearity  between  Z'  and  the  perturbations  we  may 
write  for  the  Laplace  Transform  of  the  change  in  Z'  force  due  to 
perturbations  u ' ,  w',  q<,  etc. 


2  "  °Zu  U  +  GZw^  +  GZq  q  +  —  [59J 

By  taking  the  Laplace  Transform  of  Equations  [54]  and  [55]  and 
replacing  terms  like  (z.-D-  +  Zw  '  )* '  by  w  ,  we  obtain  Equa¬ 
tion  [62]  and  Equation  [63].  If  the  motion  is  a  slowly  changing 
function  of  time,  as  is  often  the  case  in  stability  calculations, 
then  Equation  [57]  would  be  a  sufficiently  good  approximation  to 
Zi '  and  GZw  =  V  +  V  S-  Howeverh  when  an  oscillatory  mode  of 

high  frequency  exists  it  may  be  necessary  to  obtain  a  more  precise 
value  for  G^ .  The  use  of  stability  derivatives  implies  that  a 

more  accurate  approximation  to  z'  is  given  by 


Zi  '  =  +  Z  'D'+Z..'D,a+ - )w 1 

and  the  resulting  transfer  function  would  be 

g7w  =  Z  '  +  Z.  's  +  Z..  's2  +  --- 
m  w  w  v5 

Although  an  infinite  series  such  as  this  can  sometimes  give  a 
correct  value  for  G^  there  are  other  times  when  the 

not  convergent  (see  Reference  2). 


[60] 


[61] 


series  is 
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— - g.9HATI0NS  OF  MOTION  OF  AN  ELASTIC  HYDROFOIL  BOAT 

1-  Introduction 

In  the  preceding  sections  the  equations  of  motion  of 
the  hydrofoil  boat  have  been  derived  on  the  assumption  that  the 
structure  is  perfectly  rigid.  On  this  basis  the  externally 
applied  hydrodynamic  forces  are  put  in  equilibrium  with  the  hy¬ 
drodynamic  and  inertial  forces  resulting  from  the  translational 
and  rotational  motions  of  the  boat  as  a  rigid  body.  Though  this 
assumption  is  often  valid  there  is  no  assurance  that  this  is  al¬ 
ways  so,  especially  in  view  of  recent  tendencies  toward  lighter 
construction  and  higher  speeds.  Deformation  of  the  hull  and  foil- 
strut  system  may  induce  additional  hydrodynamic  forces  which 
affect  the  overall  response  of  the  boat.  In  addition,  if  the 
dynamic  response  is  a  result  of  rapidly  applied  external  forces 
such  as  those  applied  in  a  rough  sea  the  boat  will  not  only  be 
caused  to  rotate  and  translate  but  structural  vibrations  may  be 
Induced  as  well.  The  latter  may  have  a  significant  effect  upon 
the  internal  stress  distribution  of  the  structure. 

Most  of  the  important  hydroelastic  effects  on  stability 
and  control  may  be  taken  into  account  simply  by  modifying  the 
hydrodynamic  derivatives.  Thus  if  the  frequencies  of  the  rigid 
boat  motions  are  much  smaller  than  the  natural  frequencies  of 
the  elastic  deformation  modes  then  it  may  be  assumed  that  the 
hydrodynamic  loads  change  slowly  enout  so  that  the  deformation 
of  the  structure  takes  place  in  the  same  manner  as  would  be  the 
case  if  the  instantaneous  load  were  applied  as  a  static  load. 

Thus  a  change  in  load  produces  a  proportional  change  in  structure 
shape  which  in  turn  changes  the  load.  Examples  of  stability  de¬ 
rivatives  which  take  this  effect  into  account  are  given  in  Re¬ 
ferences  3-8. 
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When  the  separation  in  frequency  between  the  elastic  de¬ 
grees  of  freedom  and  the  rigid  body  motions  is  not  large  then 
significant  coupling  may  occur  between  the  two.  In  that  case  a 
dynamic  analysis  is  required.  It  was  seen  that  the  rigid  body 
equations  described  the  motion  of  the  boat  in  terms  of  six  de¬ 
grees  of  freedom:  pitch,  roll,  yaw,  surge,  sway  and  heave.  When 
the  elastic  properties  of  the  boat  are  taken  into  account  also, 
this  amounts  to  introducing  an  infinite  number  of  additional  de¬ 
grees  of  freedom.  In  practice  however,  the  elastic  deformations 
may  be  approximated  very  closely  by  a  finite  number  of  degrees 
of  freedom.  Three  methods  are  described  in  References  4  and  5 
for  reducing  a  continuous  system  to  one  with  a  finite  number  of 
degrees  of  freedom.  In  the  first,  the  deformation  is  assumed 
as  a  super-position  of  a  finite  number  of  natural  modes  of  the 
structure.  These  modes  are  obtained  from  the  eige  values  and 
eigen  functions  of  the  homogeneous  equations  of  motion  of  the  un¬ 
restrained  boat.  In  the  second,  the  deformation  is  taken  as  a 
super-position  of  a  finite  number  of  assumed  mode  shapes,  and  in 
the  third,  the  deformation  is  described  by  the  deflections  at  a 
number  of  discrete  points  on  the  surface  of  the  structure.  This 
last  method  is  sometimes  called  the  collocation  approach  and  is 
perhaps  to  be  preferred  where  automatic  computing  machines  are 
available.  This  method  is  sometimes  used  in  combination  with 
the  method  of  assumed  modes.  The  foregoing  methods  are  briefly 
outlined  below. 

2,  The  Method  of  Normal  Ixodes* 

Let  the  position  of  a  mass  element  6  of  the  boat  in 

m 

the  reference  steady  state  be  given  by  xq,  yQ,  zq  and  in  its  per¬ 
turbed  position  at  time  t  by  x,  y,  z  measured  in  the  stability 


*  - - 

This  treatment  follows  that  of  Reference  14. 
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axis  system  with  origin  at  the  mass  center.  Then  the  elastic 
displacements  of  the  elements  x-Xq,  y-yQ,  z-zq  may  be  described 
by  the  following  equations 


x-x 

0 

=  2 

1 

f 

n 

(v 

V 

z  )  e 

0 '  n 

(t) 

y-yo 

00 

=  2 

1 

Sn 

(v 

V 

z  )  e 

0  n 

(t) 

z-z 

0 

00 

=  2 

1 

h 

n 

(v 

V 

z  )  e 

0 '  n 

(t) 

where  f^  §n->  z^  give  the  shape  of  the  11-th  normal  mode  and  e 

n 

is  the  generalized  coordinate  which  gives  the  displacement  of 
that  mode.  An  important  property  of  normal  modes  is  that  in  each 
normal  mode  of  free  vibration  the  resultant  angular  and  linear 
momenta  are  zero.  Hence,  the  resulting  elastic  vibrations  de¬ 
scribed  by  the  superposition  of  normal  modes  produces  no  inertial 
coupling  to  the  Euler  equations  of  motion.  Nevertheless  the  de¬ 
formations  of  the  boat  will  produce  additional  perturbation  hy¬ 
drodynamic  forces  which  must  be  added  to  those  due  to  the  rigid 
body  motions.  Thus  the  linearized  hydrodynamic  perturbation  Z- 
force,  for  example,  would  become,  on  summing  over  n 


Z  —  Z  u+Z  w 
u  w 


+  2  Z. 
en 


e 

n 


+  2  Z. 

en 


e 

n 


+  2  Z.. 
en 


e 

n 


[65] 


Similar  additions  would  be  made  to  the  other  force  and  moment 
equations . 

Hit  additional  equations  of  equilibrium  in  the  elastic  de¬ 
grees  of  freedom  in  terms  of  the  generalized  coordinates  e  are 

n 

obtained  with  the  aid  of  Lagrange's  equations  of  motion  in  a 


movj  ng 
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frame  of  reference  (see  Reference  l) 


d  St  St  Su 

_  —  _  —  +  —  =  F 

dt  Se  Se  n 


[66] 


-n 


n 


n 


where  T  is  the  kinetic  energy  of  the  system  relative  to 
the  moving  frame  of  reference 

U  is  the  strain  energy 


F  -!* 
n  oe 

n 


W  =  work  done  on  the  system  by  the  external  forces 
(including  the  inertial  forces,  due  to  the  motion 
of  the  frame  of  reference,  which  act  on  the  system). 
Since  there  is  no  inertial  or  elastic  coupling  between  normal 
modes  the  above  equations  are  coupled  only  through  the  hydrodynamic 
contributions  to  the  F's.  Thus  the  left  hand  side  of  Equation  [66] 
for  each  normal  mode  is  Independent  of  any  of  the  other  normal 
modes.  The  kinetic  energy  of  the  n-th  normal  mode  Is  therefore 


where  I n  -  /  +  gR2  +  h^2 }  dm,  the  integral  over  the  entire 

boat,  is  the  generalized  inertia  in  the  n— th  mode.  For  a  system 
vibrating  in  an  undamped  natural  mode  the  kinetic  energy  is  a 
maximum  and  the  strain  energy  is  zero  when  the  deformation  is 
zero.  When  all  elements  are  at  the  extreme  position  the  strain 
energy  is  a  maximum  and  the  kinetic  energy  is  zero.  In  this  po¬ 
sition  the  strain  energy  is  equal  to  the  maximum  kinetic  energy. 
Thus  if  =  a  sin  co^t,  where  is  the  natural  frequency  of  the 


n-th  mode. 


then 


T 

max 


I  03  2  a 
n 


U 

max 


[68] 
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Since  the  stress-strain  relation  is  assumed  to  be  linear 

U  =  5  k  e  2,  where  k  is  a  constant,  and  U  =  Jr  k  a2 .  Combining 
n  max 

with  Equation  [68]  k  =  I  cn  2  and 

n  n 


i  I 


2  ±  03 

n  n 


e 


n 


2 


Inserting  Equations  [67]  and  [69]  into  Equation  [66]  gives 


[69] 


I  e  +  I  03  2  e 
n  n  n  n  n 


P 

n 


[70] 


In  order  to  determine  F  it  is  necessary  to  evaluate  the 

n 

work  done  on  the  system,  during  a  virtual  displacement  6e  by  all 

the  external  forces  acting  on  the  system.  Since  there  is  no  in¬ 
ertial  coupling  between  the  rigid  body  motions  and  the  normal 
elastic  modes,  the  work  done  on  each  of  these  modes  by  the  in¬ 
ertial  forces  associated  with  non-uniform  motion  of  the  frame 
of  reference  is  zero.  The  remaining  contribution  to  the  work 
done  is  that  due  to  the  hydrodynamic  forces.  Let  the  local 
normal -pres sure  perturbation  on  an  element  dS  of  the  foil-strut  ’ 
system  be  p(x  ,y  , z  )  and  let  the  local  outward  unit  normal 

vector  be  n  (n  ,n  ,n  ).  Then  the  work  done  by  the  hydrodynamic 
x  y  z 

forces  in  a  virtual  displacement  is 

5  =  -  fp  rT  •  (  ~r  -  "r o  )  dS  [  71  ] 

where  the  integral  is  over  the  whole  strut-foil  system,  and 
r  -  rQ  is  the  vector  displacement  at  dS. 

’r-"r=z(if+jg+kh)6e  [72] 

o  v  n  n  n'  n 

where  "i,  "j"  and  "k  are  unit  vectors  in  the  x,  y,  z  directions. 
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Combining  Equations  [ 71 ]  and  [72]  gives  for  P 

n 


Now  p  is  a  function  of  all  the  generalized  coordinates  (includ¬ 
ing  the  rigid  body  perturbations).  The  result  is  that  F  is  a 

n 

linear  function  of  all  these  variables  which  may  be  expressed 
in  terms  of  a  set  of  generalized  hydrodynamic  derivatives  (or 
alternatively  hydrodynamic  transfer  functions). 


dF 

F  =  ~ 
n  ou 
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[  73a  ] 


[73b] 


In  practice,  only  the  important  derivatives  would  be  retained. 
Some,  examples  of  these  derivatives  are  given  in  Reference  3.  The 
equations  may  be  non-dimensionalized  in  the  same  manner  as  the 
rigid  body  equations. 


3.  Method  of  Assumed  Modes 

When  the  elastic  deformation  is  described  by  a 
sum  of  artificially  selected  deformation  shapes  or  modes,  the 
resulting  equations  of  motion  exhibit  inertial  and  elastic  coupl¬ 
ing  between  the  modes.  In  order  to  simplify  the  illustration 
given  below  we  assume  only  the  vertical  deflections  are  important. 
Thus  the  vertical  displacement  at  a  point  x,y  of  the  structure 
may  be  written 

n 

w(x,y,t)  h1  (x,y)  e±(t) 

i=l 

where  In  describes  the  selected  mode  shapes  and  may  also  include 
the  rigid  body  modes. 
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The  kinetic  energy  of  the  system  is 


n  n 


T=|/w2dm=^^  ^ 

i=l  J=1 


m .  .  e .  e  . 

iJ  1  J 


where  the  integration  is  over  all  mass  elements  dm  of  the  struc¬ 


ture  and 


m  =  /  h.  h.  dm 
1  j 


The  strain  energy  of  the  system  is  expressed  in  terms  of  the 
stiffness  influence  coefficient  k  (x,y,  (see  Reference  4, 

page  18),  as 


U  -  2  /  w(x,y,  t)  /  k  w(^,  t),  t )  d£dq  dx  dy=^ 


where  the  stiffness  coefficient 


n  n 


I  Z  ku  v 

1=1  J=1 


k±1  =  /  h±(x,y)  /  k  h  U,T))d§dTidx  dy 

The  stiffness  coefficient  associated  with  the  rigid  body  modes 
vanish.  Inserting  T  and  U  in  Lagrange's  equations  gives 

t  L  +  t  h 3  ej  '  h  <1=1-2  — [ 

J=1  j=l 

where  Q  =  -  /  p  n  h  dx  dy 

In  contrast  to  Equation  [70]  it  is  seen  that  due  to  the  non-ortho¬ 
gonality  of  the  assumed  modes  the  inertial  and  elastic  coupling 
terms  m^  and  k±.  (l  /  j)  appear.  In  spite  of  this  added  complica¬ 
tion  it  has  often  been  shown  that  where  the  lifting  surfaces  have 
rigid  chordwise  sections  it  is  more  economical  to  choose  artificial 
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modes  because  of  the  greater  ease  of  computing  Q,  (see  pages  5&7~9 
of  Reference  4). 

4.  The  Collocation  Approach  Using  Matrices 

Longitudinal  Equations  of  Motion  -  The  treatment  given 
here  is  to  illustrate  briefly  the  general  nature  of  the  method. 

The  reader  is  referred  to  Reference  5  for  further  details.  It 
is  assumed  here  that,  for  the  longitudinal  motion,  only  the  trans¬ 
verse  (z)  deflections  are  of  hydrodynamic  importance.  In  the 
case  where  the  structure  is  divided  into  a  number  of  rigid  seg¬ 
ments,  the  torsion  is,  of  course,  of  hydrodynamic  importance  also, 
but  in  the  following  formulation  torsion  is  assumed  to  be  taken 
into  account  by  appropriate  choice  of  the  transverse  elastic  de¬ 
flection  fje^  with  respect  to  the  elastic  axis  (or  other  nodal 
line  for  twist).  Moreover,  by  virtue  of  the  fact  that  changes 
in  the  forward  velocity  perturbations  u  occur  very  slowly  com¬ 
pared  with  the  motions  of  the  elastic  perturbations  the  hydro¬ 
elastic  modifications  of  the  derivatives  in  the  X-force  equation 
may  be  calculated  on  a  quasi-steady  basis.  Therefore,  only  the 
Z-force  and  M-moment  equations  will  be  considered  in  the  follow¬ 
ing  treatment.  For  equilibrium  of  total  forces  along  the  z-axis 

f  'ZH  -  PB  V  dT  =  ° 

where  the  integration  is  over  the  entire  boat  and 

Z  =  time  and  space  dependent  hydrodynamic  Z-forces  per 
H 

unit  volume 

pR  =  mass  density  distribution  of  boat 

f  =  z-component  of  perturbation  acceleration  from  equili¬ 
brium  flight  position  of  any  point  on  boat 


[75) 


dx  =  a  volume  element  of  the  boat. 
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Por  equilibrium  of  X-forces  and  M-moments  respectively 

I  [XH  -  PB  g‘t]  dx  =  0  [76] 

1  -  PB  V  x'  [XH  -  PB  St1  z)  dT  '  °  1771 

where  x  and  z  are  the  x  and  z  coordinates  of  element  dr. 

We  next  consider  the  boat  to  be  subdivided  into  n  elements. 

The  z-component  of  the  acceleration  at  element  i  of  the  boat  is 
given  by 


f 


[Uq  -  w  +  q  xi 


[78] 


where  the  first  three  terms  on  the  right  are  due  to  rigid  body 
motion  and  the  fourth  term  is  due  to  elastic  deformation.  Hence 
the  Z-force  at  element  i  is 


Zj  =  ZH.  +  mi  [Uq  -  w  +  q  x± 


J 


[79J 


where  =  total  hydrodynamic  force  on  element  i 

irq  =  mass  of  element  i. 

JL 

By  considering  the  hydrofoil  boat  to  be  made  up  of  n  elements. 
Equation  [75]  may  be  rewritten  as 


n 


I 


i=0 


Ztj  +  m.  [Uq-w  +  qx. 
il  1  1 

1 


[80] 


The  difference  between  this  equation  and  that  for  the  rigid  body 
is  essentially  due  to  the  addition  of  the  elastic  deformation  term 

fj°^  both  explicitly  and  in  so  far  as  ZR_  depends  on  fj6^ 


and  its 
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spatial  and  time  derivatives.  The  elastic  deflection  f 


(e) 


may 


be  written  in  terms  of  the  boat  influence  coefficients  a.  .  as 


.(e) 


n 

I 

J=0 


a  .  Z  . 
ij  J 


where  a. 

ij 


if; 


(e) 


where  (f 


(e) 


=  [a.  .]  (Z.j 
ij  J 


and  [Z .)  are  the  column  matrices  of  the  elastic  de- 
J 

flection  distribution  and  total  force  distribution 


(including  inertial  effects)  given  by  Equation  [ 79] ■ 
Methods  for  determining  the  influence  coefficients  are  given  in 
References  4,  5.»  and  6. 

The  local  hydrodynamic  Z-force  column  matrix  may  be  con¬ 
veniently  broken  up  as  follows 


[8l] 


the  free  body  influence  coefficients  of  the  boat.,  re¬ 
present  the  elastic  deflection  at  point  i  due  to  a 
unit  force  at  point  j. 

Equation  [ 8 1 ]  represents  a  set  of  n  linear  differential  equations 
Written  in  matrix  form  these  become 


[82] 


(Z 


H, 


:(ZR  ]-  (A 

1  -iJ  dt2 


+  B. 


IJ 


It  +  Cij>trj(e)>  +  'V 


83] 


where  (ZR  ]  represents  the  contribution  of  the  rigid  body  motion 

ni 

on  the  hydrodynamic  Z-force  perturbation,  determined 
by  use  of  rigid  body  derivatives  and/or  with  the  aid 
of  the  square  matrix  P_. 

A.  .  =  hydrodynamic  added  mass  square  matrix  (this  matrix 
_ 1 J 

may  more  conveniently  be  included  in  P  if  desired). 
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B.  =  U  1  P 

i  J 

G.  .  =  P  D 
3-3 - 

P  =  hydrodynamic  square  matrix  by  which  the  angle  of 

^  0  \ 

attack  column  matrix  (a'  ')  must  be  multiplied  to 
obtain  the  corresponding  column  Z-force  matrix.  .P 
may  contain  the  differential  time  operator  5/Bt. 

D  is  a  streamwise  differentiation  square  matrix  defined 
in  such  a  way  that 


sil- 


} 


{a(e)]  =  (  |  1^-  +  D  )  [f(e)] 


[Zv  }  =  hydroelastic  effect  on  hydrodynamic  forces  due  to 

1 

control  surface  input  and  seaway.  These  sources  of 
Z-force  will  in  general  depend  on  the  elastic  deflec¬ 
tions  as  well  as  rigid  body  motions.  (Even  when  the 
control  system  is  assumed  known,  the  sensing  element 
of  an  autopilot  may  be  located  on  a  flexible  member. ) 


The  determination  of  the  terms  A_,  B  and  C  (the  subscript  ij  has 
been  dropped  for  convenience ),  which  represent  the  effect  of  vary¬ 
ing  local  hull  and  foil  distortion,  velocity  and  acceleration 
in  the  elastic  degrees  of  freedom, in  general  would  involve  the  so¬ 
lution  of  integral  equations.  However,  for  approximate  analysis, 
the  use* of  strip  theory  can  lead  to  useful  results  (see  References 
4-6). 


Combining  Equations  [79L  [82]  and  [83]  yields 


[f(e)]=[a 


f 


±j]  ^(ZR  ^}+m  [Uq-w+qx±]-[  C+B  ^ 


+  (m+A)^f](f(e))  +  {ZF  If  [8<l] 
”  “  St2  EJ 
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where  m  is  the  diagonal  mass  distribution  matrix. 

Equations  {84]  represent  a  set  of  n  equations  which  describe  the 
motion  of  the  boat  in  its  elastic  degrees  of  freedom.  Equation 
[ 80]  represents  the  Z-force  equilibrium  equation  of  the  entire 
boat  including  the  contribution  of  the  elastic  degrees  of  freedom. 

The  force  contribution,  from  the  elastic  degrees  of  freedom, 

Z^e\  to  the  total  Z-force  is  from  Equations  [80]  and  [83] 

Z^  =  -  1'  |  [  (m+A)  D2  +  B  D+C]  {f  +  {Z^,}  J  [85] 

where  1_'  is  a  unit  row  matrix  which  operates  so  as  to  sum 

all  the  terms  of  the  column  matrix  on  its  right. 

Since  the  hydroelastic  forces  in  the  x-direction  were  taken  as 
quasistatic  they  will  already  be  included  in  the  rigid  body 

(  0  ) 

equations  of  motion  and  since  it  is  assumed  that  is  negli- 

(  0  ) 

gible  compared  to  f}  the  contribution  of  the  elastic-degrees 
of  freedom  to  the  rigid  body  moment  are  obtained  by  multiplying 

f  0  ) 

the  terms  in  f^  '  by  -x^,  the  longitudinal  location  of  the  center 
of  pressure  of  the  elastic  element  and  summing.  Thus 

=  x'  |  [  (m+A)  d2  +b  d  +  c]  (f^)  +  {ZE)J  [86] 

where  x'  Is  a  row  matrix  of  x. 

When  Equation  [85  is  added  to  the  rigid  body  Z-force  equation 
(Equation  [51-b])  and  Equation  [86]  is  added  to  the  rigid  body  M- 
moment  equation  ([Equation  [51-c]),  the  resulting  equations,  to¬ 
gether  with  Equation  [51-a]  and  Equations  [84],  constitute  a  set 
of  ordinary  linear  differential  equations  for  the  longitudinal 
motion  of  the  elastic*"  hydrofoil  boat. 
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Lateral  Equations  of  Motion  -  The  elastic  deforma¬ 
tions  of  interest  in  connection  with  the  lateral  stability  are 
of  three  types: 

(a)  The  static,  symmetric  deflections  in  the  z-direction 
which  affect  the  equilibrium  values  of  the  angle  of  attack  and 
dihedral  distributions  of  the  hydrofoils. 

(b)  The  dynamic  lateral  deflections  in  the  y-direction 
which  affect  ail  of  the  lateral  stability  derivatives. 

(c)  The  dynamic,  anti -symmetric  deflections  in  the  z- 
direction  which  are  produced  by  the  lateral  motions  but  essen¬ 
tially  affect  only  the  rolling  moment. 

The  deflections  of  type  (a)  may  be  taken  Into  account  In  the  rigid 
body  equations  directly  and  need  not  be  considered  further  here 
The  deflections  of  type  (b)  and  (c)  will  be  treated  separately 

below  for  convenience.  The  deflections  of  type  (b )  win  be  con- 
sidered  first. 

The  lateral  force  matrix  may  be  written  in  an  analogous 
manner  to  Equations  [79]  and  [83] 

g<t>} 


'Yl]  =  {YrJ  ~  B  tv  +  Ur  -  p  z.  +  r  x± 


[  —  +  |t"  +  (m+Ai  ) 


at 


ill  ]  j-f(e)i  ,  fv  , 
2  J  +  [XElJ 


[87-1 


where  Cx  ,  $x_  ,  Ax_  etc.  have  corresponding  meaning  to  those  in 
the  longitudinal  motion  case  but  in  general  have  different  values. 

Also  (ff >)  represents  the  lateral  elastic  deflection  column  matrix. 

We  let  [aj  denote  the  structural  square  influence  matrix 
for  lateral  elastic  deflections  so  that 

=  [ax  ]  [Yx  ] 


[88] 
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Combining  [87]  and  [88]  gives  the  equations  of  motion  for  the 
lateral  elastic  deflections. 

The  contribution  of  the  effect  of  the  lateral  elastic  de¬ 
flections  to  the  Y-force j  K-moment  and  N-moment  rigid  body  equa¬ 
tions  are  obtained  from  Equation  [87]  as 

Yie)  =  -  1'  |  [(m+Ai_)  D2  +  Bi_  D  +  Ci  ]  [f;{e))  +  fy^  ]j  [89] 

Kie^  =  zx'  ([  (m+Ai  )  if  +  bx  D+Cx]  [f/e^]  +  fy  ]  1 

—  \ -  —  —  1  Ei  /  [90] 

Ni6^  =  “  V  |  [(m+Au)  D2  +  D  +  Ci_]  [f/e^]  +  [Ye  ]J  [9l] 

where  and  x^'  are  row  matrices  of  z  and  x  respectively. 

For  the  antisymmetric  elastic  effects  the  only  contribu¬ 
tion  from  the  rigid  body  motion  are  the  pure  roll  terms.  The 
force  distribution  is  given  by  the  column  matrix 

£Z2]=(zr  )-m[py.}-[C2+  D  +  (m  +  A2 )  D2 ]  [fje)}  +{Z  ]  [92] 

2 

where  (ZR  ]  represents  the  antisymmetric  vertical  component  of 

the  rigid  body  hydrodynamic  force  distribution 
>  h2  are  ttie  antisymmetric  hydrodynamic  influence 
ma  trices 

(e ) 

[f2  ]  is  the  z-component  of  the  elastic  deflection. 

We  let  [a2]  denote  the  structural  influence  matrix  for  the  anti¬ 
symmetric  elastic  deflections  so  that 
(f2(e)}  =  [a2]  (Z2) 


Combining  Equations  [92]  and  [93]  gives  the  equations  of  motion. 


[93] 
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The  total  rolling  moment  contributed  by  fje^  to  the  rigid  body 
K-moment  equation  is 


[  (m  +  A^ )  D2  +  D  +  Cg_]  [f2(e)]  +  (z^jJ 


[94] 


The  contributions  of  the  elastic  degrees  of  freedom  which  must  be 
added  to  the  rigid  body  equations  are  finally  E(e)  +  E(e) 

/v  J  i  s 

and  Nie  respectively. 


Quasi static  Equations  -  When  B  and  (m+A)  —  are 

8t2 

sufficiently  small  as  they  usually  will  be  when  the  motions  of 
the  boat  as  a  whole  are  of  primary  interest  (this  may  not  be  true 
for  high  speed  motions  in  head  seas)  then  it  is  permissible  to 
drop  these  terms  in  Equations  [84],  This  amounts  to  saying  that 
the  elastic  deflections  are  in  phase  with  the  loads  that  produce 

them.  Solving  for  (f(e)} 


} 

j 

where  I  denotes  the  unit  diagonal  matrix. 

Now  from  Equation  [85]  we  have,  on  combining  with  Equation  [95] 

[Z  }+m  (Uq-w+q  x  }  +  (Z  })  [96] 

-  1  ^EjJ  J  > 

On  adding  the  above  value  for  Z(e)  to  the  right  hand  side  of  the 
rigid  body  Z-force  equation  we  obtain  the  modified  rigid  body  Z- 
force  equation  in  which  quasi-static  hydroelastic  effects  are  in¬ 
cluded.  An  analogous  procedure  is  followed  in  modifying  the  other 
rigid  body  equations. 
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Steady  Flight  Hydroelastic  Effects  -  The  steady  flight 
hyd,roelastic  problem  Is  a  special  case  of  the  above  problem.  Ne¬ 
glecting  the  rigid  body  perturbation  terms  In  u,  v,  w,  etc.  and 
denoting  the  static  hydrodynamic  rigid  body  Z~force  distribution 
(r ) 

by  [Z  '}  and  the  weight  distribution  by  fW)  then  the  Z-force  dis¬ 
tribution  modified  by  the  elastic  deformation  of  the  foils  and 
boat  is 

(Z)  =  -  c  (f(eh  +  (z(r>) 

=  -  £  (l+[a1J]  Cf  [a.j]({Z(r))  -  CW) )  +  [Z(r)]  [ 

Similarly  static  control  effectiveness  may  be  calculated. 

If  we  put  [Zg  '}  as  the  control  rate  distribution  for  the  rigid 

control  then  the  control  rate  distribution  modified  by  the  elastic 
deformation  (neglecting  deformations  due  to  weight  distribution)  is 

(Z  )  =  -  C  (f(e))  +  (z^} 
o  6 

=  -£(1+  [a,  J  £)_1  [a  ,HZ^r))  +  (Z^r))  [ 

-J-  J  1  J  o  o 

5.  Applications 

An  attempt  has  been  made  in  the  foregoing  sections  to 
present  the  essential  features  of  the  three  principle  methods  of 
hydroelastic  analysis,  the  normal  mode,  the  assumed  mode,  and  the 
collocation  methods.  The  treatment  has  been  necessarily  brief 
and  sketchy  and  the  reader  is  referred  to  the  standard  references 
for  fuller  details  and  examples  of  their  application.  In  the 
first  ttoo  methods  the  deformation  is  assumed  as  a  super-position 
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of  a  finite  number  of  normal  or  assumed  modes  of  the  structure 
and  in  the  third  the  deformation  is  described  by  the  deflections 
at  a  number  of  discrete  points  on  the  structure.  Most  of  the 
methods  used  in  the  past  have  been  based  on  the  modal  description. 
However,  for  very  complex  configurations,  and  with  the  increasing 
speed  and  capacity  of  digital  computers  the  collocation  approach 
using  matrices  may  prove  more  feasible. 

Static  Stability  and  Control  -  Quite  obviously 

static  stability  derivatives  such  as  Z  ,  M  ,  etc.  can  be  affected 

w  w 

by  the  variation  in  the  elastic  deflection  of  hydrofoils  and 
struts  with  load.  The  importance  of  the  possibilities  of  diver¬ 
gence  of  struts  and  hydrofoils  cannot  be  overemphasized.  The 
hydroelastic  effects  on  control  effectiveness  and  reversal  is, 
of  course,  also  an  important  design  consideration.  Examples  of 
such  computations  appear  in  many  references  (see  References  4-8). 

Dynamic  Stability  and  Control  -  As  stated  earlier, 
if  the  first  natural  bending  frequency  of  the  strut  or  hydrofoil 
(taking  into  account  the  added  mass  of  the  water)  is  relatively 
large  compared  to  those  of  the  rigid  body  modes,  the  deflections 
of  the  elastic  structure  are  essentially  in  phase  with  the  per¬ 
turbation  loads  and  the  rigid  body  stability  derivatives  need  be 
corrected  only  for  the  static  elastic  deflection  produced  by  the 
loads.  If  the  first  bending  frequency  is  approximately  equal  to 
that  of  one  of  the  rigid  body  modes  then  the  correction  to  the 
rigid  body  stability  derivatives  should  include  the  dynamic  hy¬ 
droelastic  effects.  Many  examples  of  the  application  of  aero- 
elastic  methods  to  such  problems  in  the  case  of  aircraft  are  given 


HYDRONAUTICS,  Incorporated 


-56- 

in  the  aeronautical  literature  (see  Reference  5). 

Response  to  External  Forces  -  When  the  dynamic  re¬ 
sponse  of  a  hydrofoil  boat  is  a  result  of  rapidly  applied  forces 
such  as  those  occurring  in  a  rough  sea,  or  the  abrupt  deflection 
of  control  surfaces,  the  boat  is  caused  not  only  to  translate  and 
rotate  but  structural-vibrations  may  occur.  The  latter  may  have 
a  significant  effect  on  the  internal  stress  distribution  in  the 
structure  and  on  the  over-all  response  of  the  boat.  This  prob¬ 
lem  has  received  wide  attention  in  the  aircraft  field  in  relation 
to  loading  and  control  problems  in  rough  weather,  gust  alleviation, 
landing  loads,  etc.  Examples  of  applications  of  aeroelastic 
analysis  to  such  problems  are  contained  in  many  references  (see 
References  1J—  7). 


Flutter  -  Flutter  is  a  dynamic  instability  which  is 
brought  about  by  the  hydroelastic  properties  of  the  craft.  It 
may ~  involve“1rhe  unrestrained  motions  of  the  entire  craft.  This 
is  likely  to  be  the  case  when  the  first  bending  frequency  is  ap¬ 
proximately  equal  to  one  of  the  rigid  body  modes.  It  is  then 
necessary  to  include  the  rigid  body  modes  in  the  flutter  deter¬ 
mination.  When  the  inertia  of  the  strut  or  hydrofoil  is  of  the 
same  order  of  magnitude  as  the  hull,  the  rigid  body  modes  should 
be  included.  For  example  it  may  be  necessary  to  include  the  rigid 
body  roll  degree  of  freedom  in  an  analysis  of  strut  flutter.  It 
is  perhaps  much  more  likely,  however,  that  particular  modes  will 
concentrate  In  the  struts  or  hydrofoils.  When  this  occurs  one 
can  study  the  stability  of  such  modes  without  consideration  of 
the  rigid  body  modes.  One  may,  for  example,  regard  the  hydrofoil 
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support  as  stationary  and  rigid. 

It  was  long  considered  that  flutter  of  subcavitating 
^hydrofoils  was  extremely  unlikely  because  of  the  low  structural- 
to-fluid  density  ratio  (9).  (10).  More  recently,  however,  in¬ 
vestigators  (8),  (11),  (12),  (15),  and  (16)  have  demonstrated 

that  the  flutter  of  subcavitating  struts  and  hydrofoils  can  be  a 
serious  problem.  In  addition  Kaplan  and  Henry  (13)  have  shown 
that  the  same  may  also  be  true  for  supercavi tating  hydrofoils. 
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